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Measurements of the equivalent parallel conductance oj metal-insulator- 
semiconductor (MIS) capacitors are shown to give more detailed and 
accurate information about interface states than capacitance measurements. 
Experimental techniques and methods of analysis are described. From the 
results of the conductance technique, a realistic characterization of the 
Si-Si0 2 interface is developed. Salient features are: A continuum of states 
is found across the band gap of the silicon. Capture cross sections for holes 
and electrons are independent of energy over large portions of the band gap. 
The surface potential is subject to statistical fluctuations arising from 
various sources. The dominant contribution in the samples measured arises 
from a random distribution of surface charge. The fluctuating surface 
potential causes a dispersion of interface state time constants in the deple- 
tion region. In the weak inversion region the dispersion is eliminated bj 
interaction between interface states and the minority carrier band. A single 
time constant results. From the experimentally established facts, equivalent 
circuits accurately describing the measurements are constructed. 

I. INTRODUCTION 

The electrical properties of semiconductor surfaces have been 
studied for a long time. Until recently most investigators were con- 
cerned with either etched surfaces covered by a thin natural oxide 
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layer or with ultraclean surfaces that exist only in high vacuum. 
Development of improved oxidation techniques for silicon have made 
it now possible to study the electrical phenomena occurring at the 
interface between silicon and silicon dioxide. These studies resulted in 
the recognition that a semiconductor-insulator interface behaves dif- 
ferently in many ways from a "bare" surface. Of practical importance 
is the fact that an oxide-covered surface is more stable and can be 
made electrically more perfect than an unprotected surface. 

The difference between oxidized and bare surfaces will be listed here 
briefly. 

(i) There is no charge exchange with states at the air oxide inter- 
face. This type of state, the "slow surface state," has a long time 
constant and is encountered on etched surfaces. For thick oxide layers, 
only states at the oxide-semiconductor interface have to be considered. 
They are equivalent to "fast states" and will be called interface states 
in this paper. In addition to interface states, there may be traps within 
the oxide. 1 - 2 Not much is known about these traps and they will not 
be a major topic of this paper. 

(ii) All oxidized silicon surfaces contain a so-called surface charge. 
It consists of positive charges residing close to the interface. The 
surface charge originates from two sources, alkali ions 3 and built-in 
charges. 4 - 6 Alkali ions can drift through the oxide at fairly low tem- 
peratures and thus be distinguished from built-in or residual charges 
which are fixed and have a total density which depends on oxidation 
rate and crystal orientation. 6 The present paper deals with interface 
states in alkali-free systems. Reference to residual surface charge will 
be necessary in connection with their influence on the uniformity of 
surface potential. 

(Hi) Thermally oxidized Si-Si0 2 interfaces are characterized by a 
much lower density of states than "bare" surfaces. The density of 
states over most of the silicon band gap in the oxidized systems is 
usually not much greater than 10 12 crn^-eV- 1 or smaller than 
10 10 cm-2-eV -1 depending on the method of oxide preparation. For 
"bare" surfaces, density of states can be as high as 10 14 cm^-eV" 1 . 

The electrical properties of interface states are characterized by 
their density, their position in the energy gap of the silicon, and their 
capture cross section. In addition, it is necessary to know whether they 
are of the donor or acceptor type. The most widely used tool for 
investigating these interface state properties is the metal-insulator- 
semiconductor (MIS) capacitor. 

Dispersion of the capacitance can be used to obtain information 
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about the energy distribution and density of interface states as has 
been shown by Terman. 7 The capacitance technique, however, has 
severe limitations. 8 Essentially, the difficulty is that interface state 
capacitance must be extracted from measured capacitance which con- 
sists of oxide capacitance, depletion layer capacitance, and interface 
state capacitance. This difficulty does not apply to the equivalent 
parallel conductance because conductance arises solely from the 
steady-state loss due to the capture and emission of carriers by inter- 
face states and is thus, a more direct measure of these properties. 
Conductance measurements yield more accurate and reliable results, 
particularly when the density of interface states is low as in the 
thermally oxidized system because only directly measured quantities 
are used with no matching to a theoretical curve required. Both the 
capacitance and equivalent parallel conductance as functions of 
voltage and frequency contain identical information about interface 
states. Greater inaccuracies arise in extracting this information from 
the capacitance. This is illustrated in Fig. 1 which shows capacitance 
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Fig- 1 — Capacitance and equivalent parallel conductance measured at 5 kHz 
and 100 kHz on a p-type sample having an acceptor density of 2.08 X 10 16 cm" 8 
and interface state density in the 10 11 cm"*-eV" 1 range. Interface state time con- 
stants and densities are given in Fig. 14 and Fig. 15 for this sample. 
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and equivalent parallel conductance measured at 5 kHz and 100 kHz. 
The largest capacitance spread is 14 percent while the magnitude of 
the conductance peak increases by over one order of magnitude in this 
frequency range. 

In order to evaluate MIS measurements with regard to interface 
states, it is desirable to have one-dimensional current flow perpendicu- 
lar to the interface. This condition does not apply when the silicon 
underneath the field plate and beyond is inverted. Then, the capaci- 
tance and conductance are dominated by the lateral ac current flow 
that has been explained elsewhere. 10 - " Because of the positive polarity 
of surface charge, p-type silicon is normally inverted and is therefore 
most easily investigated in the depletion-accumulation region. 

In the MIS capacitor, the density of states is normally low so that 
the largest silicon band bending possible is not limited by their 
influence. Maximum possible band bending in the MIS capacitor is 
slightly less than the band gap. In this paper, interface state properties 
are described for energies in the forbidden gap within the range 2u IS 
[see Fig. 5 (a) ] only at room temperature. 

Lehovec and Slobodskoy 12 theoretically treat losses due to steady 
state capture and emission of both minority* and majority* carriers 
by a hypothetical single level interface state for the cases of zero and 
infinite recombination-generation rate through states in the silicon 
space-charge region. No direct experimental measurement and evalua- 
tion of these losses has been made until recently. 9 The theory and 
equivalent circuits developed in Refs. 12 and 13 constitute a very gen- 
eral description of all the losses and charge storage processes possible 
in an MIS structure. The gap between this general theory and what is 
actually measured in a sodium-free MIS structure is bridged in this 
work. It can be shown by conductance measurements that: 

(i) loss in the oxide measured at room temperature is negligible, 

(ii) in the samples studied, the dominant loss is by capture and 
emission of carriers by interface states rather than states in the silicon 
space charge region, and 

(Hi) the majority carrier time constant is always the dominant one. 
These three results make it possible to study the capture and emission 
properties of the interface states for each type of carrier independently. 
It is this fact which makes the MIS capacitor a powerful tool. We 
shall divide the range of surface potential from accumulation to 

♦The term minority carrier denotes those carriers which are in a minority in 
the bulk of the semiconductor crystal. The term majority carrier denotes those 
in a majority in the bulk of the crystal. 
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inversion into four regions each of which is characterized by a different 
equivalent circuit. These regions are: depletion-accumulation, mid-gap, 
weak inversion, and strong inversion. The depletion-accumulation 
region includes those values of surface potential less than u B where u B 
is shown in Fig. 5(a). Mid-gap is the region a few kT/q volts wide 
about u D . Weak inversion is between u B and 2u B . Strong inversion 
includes those values of surface potential greater than 2u B . The follow- 
ing results have been obtained using the MIS conductance technique. 

(i) A continuum of states across the bandgap of the silicon appears 
to be characteristic of the Si-Si0 2 interface. 

(it) Capture cross sections for holes and electrons are independent of 
energy over large portions of the bandgap. 

(Hi) Time constant dispersion at each bias in the depletion-accumu- 
lation region is dominated by statistical fluctuations of surface poten- 
tial. In the samples measured, these fluctuations are caused primarily 
by a random distribution of built-in residual oxide charges and charged 
interface states over the plane of the interface. 

(iv) A single time constant at each bias is observed in the region 
of weak inversion because the minority carrier capture resistance be- 
comes negligibly small and the inversion layer provides a large lateral 
conductance across the silicon surface. 

(v) Equivalent circuits are derived for the depletion-accumulation, 
mid-gap, and weak inversion regions. (See Fig. 31) The equivalent 
circuit which applies in strong inversion has been worked out pre- 
viously. 11 - 14 This circuit was proposed in Ref. 11 for the case where 
lateral ac current flow is absent and has been verified by experiment 
using n-type samples. 14 

II. SAMPLE PREPARATION 

The sample preparation processing described in this section is in- 
tended to yield oxides which remain stable only during admittance- 
voltage measurements made at room temperature or lower. The sta- 
bility requirements for device applications would be more stringent 
than this and would require somewhat different processing. No attempt 
is made to minimize built-in charge density. Both high and low inter- 
face state densities are purposely produced. 

Measurements were made on both boron and arsenic doped silicon 
oriented in the [111] and [100] directions. Samples were prepared by 
first growing a 10-/x thick epitaxial layer having a nominal resistivity 
of 1 ohm-cm on a 0.005 ohm-cm substrate. The substrate always in- 
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corporates the same dopant as the epitaxial layer. The purpose of 
using an epitaxial layer is to minimize bulk series resistance to the 
point where it can be neglected — an important consideration when 
making conductance measurements. 

The slices about 1 cm square were removed from the epitaxial 
grower and immediately placed in the oxidation apparatus described 
in Ref. 4. An oxide typically between 500 A and 700 A thick was then 
grown at 1000°C in steam at atmospheric pressure by the bias tech- 
nique of Ref. 4. 

Another group of samples was prepared by growing about the same 
thickness of oxide the same way except in dry oxygen at 1000°C 
and atmospheric pressure. Trace quantities of H 2 were removed from 
the 2 stream by first passing it through a Deoxo* unit. The 2 was 
then dried by passing it through a dry ice acetone trap. The oxide 
should be as thin as possible without being in the tunnelling range 
to make the measured admittance correspond more closely to the 
admittance of the interface. Oxidizing the sample immediately after 
the epitaxial layer is grown insures that the silicon surface has not 
had much opportunity to gather contamination from its surroundings. 
The oxide was grown by the bias technique to obtain an oxide in which 
ionic contamination is a minimum. The reason for taking pains to 
minimize ionic contamination in the oxide is that the presence of an 
ion near an interface state may perturb its energy and capture proba- 
bility although such ions do not themselves act as interface states. 3 

Immediately after oxidation, the slice was placed in an oil diffusion 
pump evaporator and a number of chrome-gold field plates were 
evaporated through a metal mask held in intimate contact with the 
oxide surface. First, a few hundred A of chromium was deposited to 
make a strong bond to the oxide. This was followed by 2000 A of 
gold to provide good electrical contact. 

The chromium and gold sources were heated to evaporation tem- 
perature in adjacent tungsten baskets. The baskets are 15 cm above 
the surface of the sample to reduce penumbra effects. 8 Since not only 
capacitance but also equivalent parallel conductance depends on the 
area of the field plate, a large field-plate diameter is used to get a 
large admittance thereby increasing detection sensitivity. However, 
making the field-plate diameter large increases the likelihood of 
encompassing defects and gross nonuniformities under it. A satisfac- 
tory diameter with these considerations in mind can be as large as 

* This unit contains a catalyst which enables H a and 2 to react at room 
temperature producing H»0. 
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1500 ii. The oxide on the substrate side was then etched off in HF. 
Cr-Au in the same proportions as in the field plate was evaporated all 
over this side to make the back contact. The slice was placed in good 
thermal contact with a massive copper block during all the evapora- 
tion steps to keep it always at room temperature. This minimizes 
sodium migration into the oxide during evaporation. Minimizing so- 
dium contamination in this step of the processing also reduces the 
likelihood of a patchy charge, distribution underneath the field plate. 

III. ORIGIN OF THE EQUIVALENT PARALLEL CONDUCTANCE 

Two experiments are described which show that the dominant 
process causing the measured equivalent parallel conductance is 
capture and emission of carriers by interface states. Series resistance 
is not important because it can be made negligible by using epitaxial 
samples or it can be measured in strong accumulation and substracted 
from the total impedance. Normally, small-signal measurements are 
made so that harmonics of the signal frequency arising from the non- 
linearity of the charge-voltage characteristic are unimportant (see 
Section 5.1). 

3.1 Loss in the Oxide 

The first experiment shows that equivalent parallel conductance 
measured at room temperature does not arise from loss in the oxide. 
Capacitance and equivalent parallel conductance for a p-type sample 
are plotted as functions of field plate bias in Fig. 2 before and after 
aging an unstable oxide in room air at 150°C. This aging produces an 
increased positive charge density in the oxide as seen by the shift of 
the capacitance curve to higher negative bias. This means that the 
electric field in the oxide resulting from the applied bias is greater for 
a given silicon surface potential after aging. However, the conductance 
curve has not changed appreciably in shape or magnitude. The con- 
ductance curve is shifted to higher negative bias by the same amount 
as the capacitance curve. This is seen by the fact that the peak of the 
conductance occurs at the same value of capacitance before and after 
aging. Therefore, the conductance is a function of silicon surface poten- 
tial and not of field in the oxide. Therefore, it is not related to loss in 
the oxide. 

3.2 Loss in the Silicon Space-Chai-ge Region 

The most important processes which can give rise to the measured 
conductance are capture and emission of carriers by interface states or 
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by states in the silicon space-charge region. It is interesting to note 
that the MIS capacitor can be used to measure the properties of either 
type of state provided one or the other dominates the loss processes.* 
Processing determines the densities and positions in the bandgap of 
levels of each type of state. In all samples made as described in Sec- 
tion II, interface state loss is dominant over the entire measurable 
range. In all these samples, no impurities other than arsenic or boron 
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Fig. 2 — Measured capacitance and equivalent parallel conductance as func- 
tions of field plate bias before and after aging at 150° C in room air. Field plate 
diameter is 3.8 X 10 _2 cm, silicon resistivity 50 fi-cm, p-type, crystal orientation 
[111], oxide thickness 1200 A, and measurement frequency 100 kHz. 

are purposely introduced into the silicon. The densities and energy 
levels of other impurities or defects inadvertently introduced into the 
silicon during fabrication of the MIS capacitor are unknown. The 
same is true for interface states. It is best therefore, to establish by 
experiment which effect dominates. This is described next. 

The magnitude of the equivalent parallel conductance measured is 
always much smaller when the oxide is grown in steam than when the 

* Sah 15 has used the MIS capacitor to study states introduced into the silicon 
by heavy gold doping. 
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oxide is grown in dry oxygen, other parameters being the same. 4 In 
addition, the oxide can be reversibly cycled from the wet to the dry 
condition by heat treatment between 200°C and 400°C in various 
ambients. To enable a quantitative comparison, it is necessary that 
the resistivity of the silicon, oxide thickness, and field plate area be the 
same for both conditions. The easiest way to accomplish this is to 
cycle the same capacitor from the wet to a drier condition. Curves 
(a) and (b) in Fig. 3 are plots of capacitance and equivalent parallel 
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Fig. 3 — (a) capacitance and (b) equivalent parallel conductance vs field plate 
bias before drying; (c) capacitance and (d) equivalent parallel conductance vs 
field plate bias of the same capacitor after heating in dry nitrogen for 17 hours 
at 350°C. Field plate diameter is 3.8 X 10 -2 cm, silicon resistivity 1 fi-cm, p-type, 
crystal orientation rill], oxide thickness 910 A grown in stream, and measure- 
ment frequency 100 kHz. 
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conductance for a steam grown oxide on p-type silicon. Curves (c) and 
(d) in this figure are measured after heating the capacitor in dry 
nitrogen (dried by passing through a liquid nitrogen trap) for 17 hours 
at 350°C. A similar set of curves for n-type is given in Fig. 4. All the 
curves in these two figures are measured at 100 kHz. Curves similar 
to these can be obtained at all frequencies in the range investigated 
(50Hz to 500 kHz) . It is important to note only that the capacitance 
in Fig. 3 from curve (<0 associated with the conductance peak of curve 
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Fig. 4 — ( a ) capacitance and (b) equivalent parallel conductance measured vs 
field plate bias before drying; (c) capacitance and (d) equivalent parallel con- 
ductance vs field plate bias of the same capacitor after heating in dry nitrogen for 
17 hours at 350°C. Field plate diameter is 3.8 X 10- 2 cm, silicon resistivity 0.75 fi-cm, 
n-type, crystal orientation [111], oxide thickness 600 A grown in stream, and 
measurement frequency 100 kHz. 
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(d) is nearly the same as the capacitance from curve (a) associated 
with the conductance peak of curve (b). The same is true in Fig. 4. 
Therefore, all the relevant parameters before and after drying are the 
same except the magnitudes of the conductance peaks. The peak con- 
ductance of curve (b) in Fig. 3 is about one-third the peak conductance 
of curve (d). This result cannot be due to sodium contamination 
which might have been introduced into the oxide during the drying 
cycle, the initial wet oxide being sodium free (see Section II). The 
reasons for this are: (i) Sodium at the interface always remains 
charged in the range between accumulation and strong inversion. 3 
Therefore, it does not become an interface state which could influence 
the conductance, (ii) A nonuniform charge distribution in the oxide 
caused by the sodium would spread the conductance curve over a 
larger range of bias which would decrease rather than increase the 
magnitude of the conductance peak if this effect were dominant. 
(iii) The p-type curves in Fig. 3 are shifted towards more negative 
bias while the n-type curves in Fig. 4 are shifted towards less negative 
bias with drying. Introduction of sodium during drying would shift 
both curves towards more negative bias. The increase of the conduct- 
ance peak in both p and n-type is then a true drying effect. 

The fact that the oxide can be reversibly cycled from the wet to 
the dry state at low temperatures rules out the possibility that a large 
density of structural defects or impurities which act as bulk states is 
introduced into the silicon during one method of oxide preparation and 
not the other. All of this strongly suggests that the measured equiv- 
alent parallel conductance is due mainly to losses in states which 
reside on the oxide side of the interface and exchange charge at 50 Hz 
to 500 kHz with the silicon because only the properties of such states 
could be so markedly altered by low temperature wet-dry cycling. 

The most convincing evidence, however, is the dependence of inter- 
face state time constants on surface potential to be discussed later in 
this paper. 

3.3 Loss Due to Minority Carrier Transitions 

Minority carrier density can respond to the applied signal frequency 
by diffusion from the bulk to the interface and by recombination- 
generation processes through interface states and states in the silicon 
space charge region. These processes in the samples measured are too 
slow to allow the minority carrier density to follow 50 Hz to 500 
kHz. 11 - 14 It was indicated in the preceding part of this section that in 
this frequency range, interface state loss is solely responsible for the 
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measured equivalent parallel conductance. In the depletion region 
where minority carrier density is orders of magnitude less than 
majority carrier density, loss due to minority carrier transitions will be 
negligible compared to loss due to majority carrier transitions. Thus, in 
this region, the conductance measured arises simply from capture and 
emission of majority carriers by those interface states within a few 
kT/q of the Fermi level. In the region of weak inversion where 
minority carrier density is orders of magnitude greater than majority 
carrier density, minority carrier time constant becomes shorter than 
majority carrier time constant. Minority carriers will follow the signal 
frequency. However, current can flow and produce a loss only by 
transitions through the interface states from the minority carrier band 
to the majority carrier band. Therefore, this process will be controlled 
by the majority carrier time constant. This case will be treated in 
detail in Section IV. The important point is that the loss is determined 
entirely by the majority carrier time constant in both the depletion 
and weak inversion regions. 

IV. THEORY 

4.1 Phenomenological Description of Loss 

A continuum of interface states over the bandgap of the silicon 
appears to be characteristic of the Si-Si0 2 interface. Fig. 5(a) is a 
schematic of an MIS capacitor showing the band bending in depletion 
for n-type silicon with an arbitrary distribution of interface states in 
the silicon bandgap. Initially, let the interface states be in equilibrium 
with the silicon. Then let a small ac signal between 50 Hz and 500 kHz 
be applied. Consider the first half of the cycle which moves the con- 
duction band towards the Fermi level. This immediately increases the 
average energy of the electrons in the silicon. Because a conductance is 
observed, it is evident that the interface states do not respond im- 
mediately but lag behind. Therefore, electrons at a higher average 
energy in the silicon will be captured by interface states at a lower 
average energy. This results in an energy loss. On the other half of 
the cycle, the signal moves the valence band towards the Fermi level. 
Electrons in filled interface states now will be at a higher average 
energy than electrons in the silicon. As the electrons are emitted by the 
interface states into the silicon they will lose energy again. Thus, there 
will be an energy loss on both halves of the cycle which must be sup- 
plied by the signal source. The energy required for transitions between 
the band edge and an interface state is much higher than this and is 
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Fig. 5 — (a) Schematic showing band bending in depletion at silicon surface in 
an n-type MIS capacitor. An arbitrary distribution of interface states is shown 
for illustrative purposes. Majority carrier transitions to two levels within kT/q 
of the Fermi level are also shown. /„ is the Fermi function, (b) Band bending 
in depletion showing point at which Fermi level crosses a shallow impurity level. 
This impurity level is hypothetical. W is space charge width, (c) Band bending 
in weak inversion showing point at which Fermi level crosses a hypothetical im- 
purity level located at mid-gap. 



supplied by phonons. Only majority carrier capture and emission is 
considered here because this is the dominant loss mechanism for 
reasons given in Section III. 

After a disturbance has been applied, the interface states will reach 
equilibrium with the silicon exponentially with a characteristic time, 
t. The interface states are also capable of storing charge. Thus, there 
will be a capacitance, C„, associated with them which is proportional to 
their density. The conductance related to the loss is G 8 = C«/t. It 
should be noted that G 8 is independent of signal amplitude as long as 
C 8 and t are independent of it. This will be the case for small-signal 
amplitude. Interface state t and C s can then be extracted directly from 
a measurement of equivalent parallel conductance. Capture cross 
section is extracted from t and density of states from C s as will be 
shown in Section VII. 
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4.2 Admittance of a Single Level State 

Our purpose is to derive expressions for the admittance of the MIS 
capacitor as a function of bias and frequency, and to obtain the cor- 
responding equivalent circuits which fit in detail to experimental obser- 
vations. The capacitance-voltage characteristics for the MIS capacitor 
without interface states has been described. 10 - 10 - 17 Capacitance dis- 
persion caused by interface state loss was measured by Terman 7 to 
obtain interface state densities and time constants. Lehovec et al, 12 
Bertz, 18 and Sah 15 developed equivalent circuits for a hypothetical 
single-level interface state from Shockley-Read statistics. 18 The equiv- 
alent circuits given in this previous work are either too general 12 - 13 - 15 
or the measurement technique too inaccurate 7 to give a real and com- 
plete description of interface state properties. To describe the MIS 
conductance technique and interpret the measurements in the depletion 
and weak inversion regions, we shall rederive the theory for the ad- 
mittance of the MIS capacitor. The starting point will be the theory 
for majority carrier capture and emission by a single-level interface 
state first calculated by Lehovec and Slobodskoy. 12 

Experimental evidence shows that only capture and emission of 
majority carriers are important when measuring in the depletion 
region. The quasi Fermi level and Fermi level are identical for major- 
ity carriers. Thus, application of an ac signal simply results in a time 
varying Fermi function. From Shockley and Read, 18 the capture rate 
of electrons, taken as majority carriers, by a single-level interface state 
is 

R n (t) = N,c„[l - Ki)]n.(t), (1) 

and the emission rate is 

G n (t) = N,e n j(t), (2) 

where N, is the density of states cm -2 ; c„ electron capture probability,* 
cm 3 /sec, e n electron emission constant, sec -1 , j(t) the Fermi function 
at time t, and n.(t) electron density at the silicon surface at time t, cm -3 . 
Net current density flowing is 

i,(t) = qN,c n [l - 1(t)]n,(t) - qN.e n 1(t), (3) 

where q is electronic charge in coulombs. Making the small-signal 
approximation, the admittance is, as is shown in Appendix A, 

* Thia capture probability is the average over all states in the conduction band 
near the silicon surface. 
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Y - = >»kf(l+M /c n n.y (4) 

where j = V — 1, u the angular frequency of the ac signal, sec -1 , 
A; Boltzman's constant, eV X coul/°K, Tthe absolute temperature, °K, 
/„ the Fermi function established by the bias, j = [1 + exp (u — u.)]' 1 , 
and n to the electron density at the silicon surface established by the 
bias, cm -3 . 

Equation (4) is the admittance of a series RC network with capac- 
itance C, = q 2 N,j (l — j„)/kT and time constant r = j„/c n n, . Separat- 
ing (4) into its real and imaginary parts, the equivalent parallel ca- 
pacitance is 

C - i-JHn , (5) 

1 + CO T 

and the equivalent parallel conductance is 

/-» C,co T ra\ 

G„ = ,,22 - (6) 

1 fur 

4.3 Equivalent Circuit of MIS Capacitor 

Let Q T be the total charge density at a given bias. Then 

Qt = Q.c + Q. + Qr, (7) 

where Q, c is the silicon space-charge density, coul/cm 2 , Q, the interface 
state charge density, coul/cm 2 , and Q f the fixed-charge density in the 
oxide, coul/cm 2 . The ac current density, i T (t), obtained by differentiating 
(7) with respect to time is 

ir(t) = i. c {t) + t.(0. (8) 

where i, c (t) and i,(t) are the ac current densities charging the silicon 
space-charge layer and the interface states, respectively, To obtain 
i, c (t) we have 

»•••«> - @r)(f ) - 

where \p,{t) = \j/, + H. is the silicon band bending in volts at time, 
t, \p, the silicon band bending established by the bias, and h\p. — 
a exp (jut). From this, dfjdt = ja d\p. . Also, dQ.Jdf, = C D the deple- 
tion layer capacitance per cm 2 . Substituting these into (9) we get 

i, t = jccCn ty. . (10) 
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It is shown in Appendix A that i,(t) given by (3) can be written 

i.(t) = Y. W. , (11) 

where Y. is denned by (4). Substituting (10) and (11) into (8), we have 

ir(t) = (jcoC D + Y.) 8+. . (12) 

From (12), it is seen that C D appears in parallel with the series RC 
network of the interface states. The voltage, v a = v„ + 8v a , applied 
to the capacitor divides between the silicon and the oxide film so that 



v„(t) =v + 8v a = fcfl) + §*- , 



(13) 



where v„ is the dc bias, 8v a = b exp (jut), and C BX is the oxide capacitance 
per cm 2 . Differentiating (13) with respect to t to get the ac terms only 



]u ov a = jco by, -+- -^ — 



Substituting for 8\[/. from (12) 



8v a = i T (t)(Z. + - 



1 



ju>C c 



(14) 



(15) 



where 



Z. = (io,C D + Y.y\ 

The bracketed term in (15) is the impedance of a circuit in which 
C ox is in series with Z, . The equivalent circuit for the MIS capacitor 
from (4), (12), and (15) is given in Fig. 6. 



c D : 



i 



Fig. 6 — Equivalent circuit of MIS capacitor for capture and emission of ma- 
jority carriers by a single level interface state. Co* is the oxide capacitance per- 
unit area, Co the depletion layer capacitance per unit area, C, the interface state 
capacitance per unit area, and R. the majority carrier capture resistance. Bulk 
generation-recombination is zero. 
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The principle of the MIS conductance technique is easily illustrated 
by this simplified equivalent circuit. The admittance of the capacitor 
is measured by a bridge across the terminals x-x. The oxide capacitance 
is measured in the region of strong accumulation. The admittance of 
the network is then converted into an impedance. The reactance of 
the oxide capacitance is subtracted from this impedance and the re- 
sulting impedance converted back into an admittance. This leaves C D 
in parallel with the series RC network of the interface states. The 
capacitance from (5) and (12) is 

C p = C n + ,_?',, , (16) 

1 + CO T 

and the equivalent parallel conductance divided by co from (6) is 

w 1 + to T 

G p is divided by w to make (17) symmetrical in <ot. It should be noted 
that (16) and (17) are equivalent to the Debye equations, 19 - 20 for a 
dielectric in which the polarization can relax exponentially with char- 
acteristic time, t, after a change in applied electric field occurs. 

Equation (16) describes the capacitance dispersion and is the basis 
of Terman's method. To extract C 8 and t from C p using (16), C D 
must be known. C D can be calculated 16 using an estimated doping 
density. The doping density is not accurately known near the silicon 
surface because of pile up or depletion of the dopant during oxidation. 
Equation (17) on the other hand depends only on the interface state 
branch of the equivalent circuit. G p /«> goes through a maximum when 
wt = 1 which gives t directly. The value of G p /o> at the maximum is 
C s /2. Thus, equivalent 'parallel conductance corrected for C ox gives C 8 
and t directly from the measured conductance. 

4.4 Admittance of an Interface State Continuum 

The interface states are observed to be comprised of many levels so 
closely spaced in energy that they cannot be distinguished as separate 
levels. Thus, they appear as a continuum over the bandgap of the 
silicon. These observations are documented in Section VII and dis- 
cussed in Section VIII and agreed with those of other investigators. 7 ' 21 
A continuum of interface states appears to be characteristic of the 
Ri-Si0 2 interface. The admittance given by (4) for a single-level state 
and the equivalent circuit in Fig. 6 must therefore be modified. 

For a continuum of states at a finite absolute temperature, capture 
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and emission of majority carriers can occur by states located within a 
few kT/q on either side of the Fermi level. This results in a time 
constant dispersion. The admittance of the continuum first calculated 
by Lehovec 22 is obtained by integrating (4) over the bandgap 

v . (q 2 \ f Nj q - ;.) d± (vx\ 

Y - = nw J (i + M /c n n. ) (18) 

Here N a8 is the density of interface states, cm- 2 eV _1 , and \p is energy 
in ev. The integrand of (18) is sharply peaked about the Fermi level 
with a width of about kT/q. Thus, (18) can be easily integrated if both 
N 88 and c„ the capture probability* do not vary very much with $ over 
a range kT/q. Experimentally, it is found (see Section VII) that 
neither vary appreciably over kT/q. Making the substitution f {l—f ) 
= {kT/q) (dfo/<fy) transforms (18) into an integral over /„. Integrat- 
ing from zero to unity yields. 

Y„ = 3E±i ln (1 + w 2 t 2 ^ + j q Kn arc tan (cor m ) (19a) 

where m 

r m = l/c»n so . (19b) 

Concentrating on the real part of (19a), we have for the continuum 
& = jJr"ln(l+«'Ti). (20) 

CO *0)T m 

At no value of bias does (20) even closely fit experimentally measured 
Gp/<a vs w curves. As will be shown in the next section, (20) is just the 
first step in deriving an expression for G p /<o which fits the experiment 
in depletion. To proceed further, it is necessary to consider separately 
the depletion region and the weak inversion region. In depletion, the 
time constant dispersion observed is far broader than given by (20). 
In weak inversion, the time constant dispersion disappears and the 
experimental G p /<* vs w curves are characterized by a single time con- 
stant at each bias. 

4.5 Depletion Region 

The model suggested here for the increased broadening of the time 
constant dispersion observed experimentally (see Section VII) in the 
depletion region assumes statistical fluctuations of surface potential 
in the plane of the interface. Majority carrier density at the silicon 



* This capture probability is now the average over all the interface states in 
the range of the integral and all states in the conduction band near the silicon 
surface. 
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surface is related to surface potential by 

, n, = N D exp u, = ?i, exp (w, — u B ) for n-type ,„... 

p.o = N A exp (— u,) = n, exp —(u, — u B ) for p-type, 

where n+ is intrinsic carrier density, enr 8 , u B = In rii/N D for n-type 
and u B = In N A /n t for p-type. u-b is the potential difference between 
raid-gap* and Fermi level in the quasi neutral region of the silicon in 
units of kT/q, N D is the ionized donor and N A the ionized acceptor 
density in the silicon, cm -3 , u 8 is the silicon band bending or surface 
potential in units of kT/q. u a — u B is the potential difference between 
mid-gap and Fermi level at the silicon surface. 

The relation between time constant and surface potential is from 
(19b) and (21). 

r m = exp — (u, — u B ) for n-type 

and CnUi (22) 

r m = exp (u, — u B ) for p-type. 

Equation (22) shows that small fluctuations in u 8 will cause large 
fluctuations in r m and thus broadened time constant dispersion. If we 
assume that the built-in charges and charged interface states are ran- 
domly distributed in the plane of the interface, the electric field at the 
silicon surface will fluctuate over the plane of the interface. This is 
shown schematically in Fig. 7. Fluctuations in electric field will cause 
corresponding fluctuations in surface potential. To express this quanti- 
tatively^ we conceptually divide the plane of the interface into a 
number of squares of equal area. The area of each square is the largest 
area within which the surface potential is uniform. We shall call this 
particular area a characteristic area. Now, the admittance of the 
continuum given by (19) can be regarded as the result of integrating 
the admittance of a single level (4) over all the levels located within 
a characteristic area rather than within the entire area of the field 
plate. The total admittance is obtained by integrating the contri- 
bution from each characteristic area over all the characteristic areas 
within the area of the field plate. This must be done for each bias and 
frequency to generate the complete family of curves. The calculated 

* Although the intrinsic Fermi level is slightly above mid-gap in silicon, this 
is neglected and intrinsic Fermi level is taken to be at mid-gap. 

t This treatment is similar to the one used by W. Shockley in Ref. 23 for 
calculations of the influence of fluctuations of doping on avalanche breakdown. 
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G p /oi curves will be compared with the measured (? p /w curves to see 
whether the model agrees with observation. For simplicity, we shall 
derive in detail only an expression for the real part of the interface 
state admittance. The steps in deriving the imaginary part are similar. 
Let P{N) be the probability that there are N built-in charges and 
charged interface states in a characteristic area. The number of 
characteristic areas which contain between N and N+dN surface 
charges is 

dp = P(N) dN. (23) 



RANDOMLY 

DISTRIBUTED 

CHARGES 



SURFACE 
FIELD-- 
LINES 



METAL FIELD 
•'PLATE 




fflffl ffl fflffifflffl QQ <3<5X3> 



-SL0 2 



DEPLETION 
LAYER 



NEUTRAL 

h SILICON 

p-TYPE 



Fig. 7 — Cross section of MIS structure illustrating the random distribution of 
built-in charges and charged interface states over plane of the interface. Field 
lines are all terminated at edge of depletion layer rather than on ionized accep- 
tor impurities in the layer for simplicity. 

The total G p /<» is G p /a> for the continuum of levels located in a charac- 
teristic area (20) times the number of characteristic areas containing 
between N and N+dN surface charges (23) integrated over all the 
characteristic areas. To integrate this product over all the characteris- 
tic areas, (20) is expressed as a function of u B using (22) and P{N) 
transformed to P(w 8 ), the probability that the surface potential in a 
characteristic area is u 8 . This transformation will be performed in two 
steps. The first transformation is 



P(Q) = P(N) dN/dQ , 



(24)< 



* See Ref. 24 for a discussion of transforming the probability density as a func- 
tion of one random variable to another. 
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where Q is the density of built-in charges and charged interface states, 
coul/cm 2 . The second transformation is 

P(u.) = P(Q) dQ/du. . (25) 

The number of characteristic areas in which the surface potential 
is between u„ and u 8 4- du„ is P{u s )du g . Multiplying this by (20) and 
integrating, the total G p /m is 

(?,/« = &* f °° — In (1 + u*rl)PM du. . (26) 

It is assumed in (26) that both interface state density and capture 
probability are uniform over several kT of bandgap energy. That these 
are reasonable assumptions can be seen from Fig. 14 and 15. 

Substituting (22) for p-type* into (26) to get (26) as a function 
of u„ 



Gv/u = iqN. t / cxp (u B — u, — u n ) 

•/ — no 



•In (1 + exp [2(w + u. - u B )]\P(u,) du, , (27) 

where u = In co/c„n, . 

The next task is to obtain P(u.) from P(N). Let N be the mean 
number of built-in charges and charged interface states in a char- 
acteristic area. When N is large, P(N) is given by the Gaussian ap- 
proximation of a Poisson distribution 

P(N) = (2iriV)"' exp [-(JV - Nf/2N]. (28) 

The characteristic area a is the ratio between the mean number of 
surface charges N in a and their mean density n: a = N/n. From this, 
the relation between N and Q required for the first transformation is 

N = aQ/q. (29) 

Combining (24), (28), and (29), we get 

P(Q) = {2™Q/q)-H«/q) exp [-«(Q - Q) 2 /2qQ], (30) 

where Q is the mean of Q. 

The transformation of P(Q) to P(u.) consists of transforming an area 
in the space made up of all the points Q into an area in the space made 
up of all the points u, . This requires a single-valued relation between 

* Thia part of the derivation is carried out for p-type because the sign of the 
variables for n-type is negative making the analysis harder to follow. 
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Q and u, . A relation between Q and u, which is not single valued is 
obtained from (7) and (13). Retaining only the dc terms in (13) and 
expressing surface potential in units of kT/q = (T , we get 

Q = Q. + Q f = C 0X (v o + u./fi) + Q, c (u.). (31)* 

Q, e (u.) is the space-charge density which from Ref. 16 is 

Q. c (u.) = (2g €li iVV0)*[exp (-**.) + «. - 1]*, (32) 

where c., is the dielectric permittivity of silicon, farads/ cm and N A 
the ionized acceptor density, cm -3 . Equation (31) is not a single-valued 
relationship between Q and u, because of (32). The difficulty in making 
the transformation from Q space to u, space imposed by (31) and (32) 
is avoided by restricting the problem to the case where the fluctuations 
Q — Q are very small. For small fluctuations, (31) can be differentiated 
assuming N A and oxide thickness to be uniform to get 

dQ = Gjz du, + dQ sc . (33) 

dQ, c is eliminated from (33) by using the relation dQ, e = (C D //3) du, = 
(e.i/TTjfl) du, where W is the space-charge width. Because Q is given 
by (31) when u, = u, , dQ can be evaluated about Q at each bias by 
setting W — W(u,). Doing this, (33) becomes 

dQ = ^^ [W{u,)C 0I + e„.] du, . (34a) 

Equation (34a) is the transformation equation we seek from an element 
of area in Q space to the corresponding element in u, space. Replacing 
dQ and du, in (34a) by the small fluctuations Q — Q and u, — u, , re- 
spectively, we get 

Q-Q = ppj^ [W(u,)C BI + e, { ](u, - u,). (34b) 

Combining (25), (30), (34a), and (34b), we get 

P(u.) = (27r<r;r 5 exp [-(«. - u.f/2<T 2 .], (35) 

where the standard deviation a, is 

_ W(*.)fi (qQV ™ 

a '~ [W(u,)C 0I +€.,] U I ■ m) 



* In this paper, Qt + Q. will be calculated from flat baud voltage shift 
neglecting work function differences. Work function differences are neglected 
because flat band voltage shifts of several volts independent of the field plate 
metal are found in the samples used. 
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Q can be calculated from (31) and (32) at u, = u, . In the samples 
measured, the total number of built-in charges was more than an order 
of magnitude greater than the total number of charged interface states. 
Therefore, the bias dependence and the variation during each cycle of 
the signal of the number of charged interface states is negligible. 

Another cause of surface potential fluctuations which shall be taken 
into account is the random distribution of ionized acceptors in the space 
charge region. This is calculated in Appendix B using a model which 
has previously been used to determine the fluctuations of breakdown 
voltage. 23,25 Essentially, the procedure is similar to the one for the 
surface charges except a characteristic cube having a side equal to W 
is used instead of a characteristic area and Q rather than N A is assumed 
to be uniformly distributed. The standard deviation from Appendix B, 
paragraph B.l is 

_ qfi[N A W(U.)'\ i [l ~ exp (-fl.)] „ 

2[W(u.)C ox + e.J ' {6 ' 

where N \ is the mean ionized acceptor density. 

A theorem in statistics states 26 that if two or more independent 
variables are normally distributed their sum is also normally distributed 
with a mean and variance which is the sum of the mean and variance 
of each variable. Using this theorem, P(u.), which includes the influence 
of both the randomly distributed surface charges and ionized acceptors 
on the fluctuations of surface potential is 

P(u.) = [2r(a] + 4)]"* exp [-(«. - *.?/2{<r\ + <£)]. (38) 

Substituting (38) into (27) 

GJ<a =hqN..[2v(<r* + •»]-* 

■f exp [-(g + ?/)] In (1 +e 2v )du, , (39) 

J -00 

where y = In coT m = u + u, — u B and z = (u, — u,) 2 /2(o 2 , + o-J). 
By similar arguments using (19), the equivalent parallel capacitance is 

C p = C D {u.) + qN. ,[2ir(a] + «£)]"* 

/ ex P [ — ( 2 + y)] arc tan (e") du, . (40) 
-00 

Depletion layer capacitance C D now becomes a function of u, . The 
function C D (u.) has the same form as in Ref. 16 as long as the fluctua- 
tions Q — Q are small. 
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It will be shown in Section VIII that the influence of the randomly 
distributed ionized acceptors on the fluctuations of surface potential 
are small in the samples measured compared to the influence of built-in 
charges. 

In summarizing the theory for the depletion region, it is shown how 
the interface state branch of the equivalent circuit is developed. The 
interface state branch of the equivalent circuit of Fig. 6 applies to a 
single level in the continuum. In each characteristic area where the 
surface potential is uniform, each level contributes a series RC network. 
Integrating the admittance of each level given by (4) over all the levels 
located in a characteristic area gives (19) and (20). Integrating again 
but this time the contribution to the total admittance from each char- 
acteristic area over all the characteristic areas gives (39) and (40). 
The interface state branch of the equivalent circuit therefore, can be 
represented by an infinite number of series RC networks connected in 
parallel as illustrated in Fig. 8(a). 

4.6 Weak Inversion 

For n-type, the silicon surface will be in weak inversion when the 
Fermi level is more than a few kT/q below mid-gap in the lower half 
of the bandgap. The minority carrier density at the silicon surface will 
now be orders of magnitude greater than the majority carrier density 

I 




Cp(w)^ 



•G P (w) 



(a) 



(b) 



Fig. 8 — (a) Equivalent circuit for depletion region showing time constant 
dispersion caused primarily by statistical fluctuations of surface potential. C » 
is the oxide capacitance per unit area, and C D the depletion layer capacitance per 
unit area. Each subnetwork consisting of C. and R. in series represents time con- 
stant of the continuum of interface states in a characteristic area. Bulk generation- 
recombination is zero, (b) Simplified version of (a) useful when extracting C B * 
from the measured admittance. C P (u) is the capacitance per unit area at a given 
bias and frequency of the distributed network in parallel with C D . G p (w) is the 
equivalent parallel conductance per unit area. 
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at the silicon surface. The minority carriers can no longer be ignored 
as in depletion and the equivalent circuit of Fig. 8(a) does not apply. 
To derive expressions for the interface state admittance we start again 
by considering a single level state. The capacitance associated with 
this state is C. = q 2 N ,j {\ — j )/kT and its time constant for majority 
carrier transitions (electrons) is t„,, = j /c„n lo as derived previously 
in paragraph 4.2. The majority carrier capture resistance from this is 

R n .. = T.../C. = (kT/q 2 )[c n n. B N.(l - /.)] _1 . (41) 

By similar reasoning, it can be shown that the minority carrier 
(holes) capture resistance is 

«,.. = (kT/q 2 )[c vV . N.UY\ (42) 

where c v is the capture probability for holes cm 3 /sec and p lo is the 
hole density cm -3 at the silicon surface when the Fermi level is at the 
interface state level. 

The equivalent circuit for a single level is given in Fig. 9(a). From 
(41) and (42) 

R«,./Rp,. = (c p /c n ) exp (u F - 3m). (43) 

Equation (43) has been obtained by substituting p, = n, exp ( — u), 
n,„ = n { exp (u), and /„ = [1 + exp (u — Up)]' 1 into (41) and (42). 
u is the potential difference between mid-gap and the interface state 
level in units of kT/q and u F is the potential difference between mid-gap 
and the Fermi level at the silicon surface in units of kT/q. There will 
be no detectable loss due to transitions of electrons and holes to and 
from interface states unless u tt u F because / (1 — /„) is so sharply 
peaked around the Fermi level. Equation (43) becomes 

Rn,./R P .. = (c P /c n ) exp ( - 2u F ) . (44) 

It is seen from (44) that R„, g will rapidly become orders of magnitude 
greater than R p% , when the Fermi level gets more than a few kT/q 
below mid-gap even if c p is one or two orders of magnitude smaller 
than c„ . Time constant for majority carrier transitions is r n ,, = C.R n ,, , 
and for minority carrier transitions it is t,,, = C,R P ,, . Thus, t„., will 
be much greater than t p ,. in weak inversion. If the period of the applied 
signal frequency is comparable to t„,„ as it is in these experiments, 
then there will be virtually no loss associated with minority carrier 
transitions. Each interface state near the Fermi level is charged and 
discharged instantaneously by minority carriers so that R Pi , tt 0. 

Because R p _„ is negligibly small, it can be replaced by a short circuit 
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R n,s^ < R P,s R n, 
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T , T T 
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Fig. 9 — (a) Equivalent circuit for a single level interface state in region of 
weak inversion. C. is the interface state capacitance per unit area, R n ,. the cap- 
ture resistance for electrons, and R P ,. the capture resistance for holes. CB repre- 
sents the conduction band and VB the valence band, (b) Equivalent circuit for 
a single level interface state for the case in which minority carriers are holes 
which respond instantaneously to the applied signal. For this representation, it 
is also necessary that the lateral conductance due to minority carriers along the 
surface is large compared to Rn,.' 1 . (c) Circuit (a) generalized to include the fact 
that there is a continuum of states across the bandgap of silicon. C/ is the in- 
version layer capacitance per unit area, (d) Circuit (c) simplified by fast response 
of minority carriers. Ci is negligible in the region of weak inversion, (e) Simplifica- 
tion of circuit (d) . 



which connects the interface state capacitance directly to the minority 
carrier (valence) band as shown in Fig. 9(b). Those interface states 
in the continuum which are located in a given characteristic area each 
have a different time constant. These interface states must be repre- 
sented by a distributed network such as Fig. 9(c). Only those levels 
in a narrow range of energy about a kT wide centered about the Fermi 
level participate in the capture and emission process. Therefore, the 
corresponding range of time constants is narrow. The large minority 
carrier density at the surface then means that the minority carrier 
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capture resistance R Vi . of every one of these levels is negligibly small. 
Thus, the capacitance corresponding to each of these levels can be 
connected directly to the minority carrier band. The large inversion 
layer conductance shorts all the characteristic areas together. The net 
result is represented in Fig. 9(d) which can be further simplified to 
Fig. 9(e). G p /o) vs to now will be characterized by a single time con- 
stant at each bias. The expression for G P /u is (17) and for C P it is (16). 
The value of C. in (16) and (17) is the sum of all the individual interface 
state capacitances in parallel as shown in Fig. 9(c). The inversion 
layer capacitance C T calculated in weak inversion using the relation 
given in Ref. 12 is negligible compared to 2C, obtained from meas- 
urement. 

If interface state density N,, is slowly varying with the position of 
the Fermi level, the statistical fluctuations of surface potential will 
have negligible effect on the total C, . However, in analyzing the meas- 
urements on n-type samples in weak inversion using (17), it is found 
that C. varies quite rapidly with position of the Fermi level (see Fig. 27). 
This implies that N, t is also rapidly varying so that the statistical 
fluctuations of surface potential must be taken into account when 
calculating the relation between N„ and total C, . Total capture con- 
ductance is the sum of the capture conductances R~) of the individual 
interface states in parallel as shown in Fig. 9(e). 

The statistical fluctuations of surface potential must be taken into 
account in calculating the total capture conductance because capture 
conductance depends on n, as well as the rapidly varying N„ . If 
desired, C. can be calculated by integrating C, for a single level over 
bandgap energy and then over surface potential using the statistical 
model developed earlier. Similarly, total capture conductance can be 
calculated by integrating capture conductance for a single level ob- 
tained from (41) over bandgap energy and then over surface potential. 
These integrals cannot be evaluated without knowing the functional 
dependence of N ,, on the position of the Fermi level. This dependence 
can be found by numerically fitting the calculations to the measure- 
ments. Although N,, and majority carrier capture probability can be 
obtained in this way from measured C, and t„ , such an analysis will 
not be carried out here. 

A comparison between the depletion and weak inversion regions 
will illustrate the processes occurring. In depletion except very near 
mid-gap, majority carrier density at the silicon surface is several 
orders of magnitude greater than minority carrier density at the silicon 
surface. Therefore, the time constant for minority carrier transitions 
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is much greater than the time constant for majority carrier transitions. 
The time constant for majority carrier transitions is comparable to 
the period of the applied signal frequency. Therefore, minority carrier 
density cannot follow the signal at all and ac current flows simply by 
transitions between interface states near the Fermi level and the 
majority carrier band. Thus, the interface state branch of the equivalent 
circuit in Fig. 8(a) can be obtained by open circuiting R v ,. in Fig. 9(c). 

In the weak inversion region except very near mid-gap, again the 
time constant for majority carrier transitions is comparable to the 
period of the signal frequency and several orders of magnitude greater 
than the time constant for minority carrier transitions. Now both 
majority and minority carrier densities respond to the signal so that 
ac current flows by generation and recombination through interface 
states near the Fermi level. This current flow is controlled by the 
majority carrier capture resistance as it is the largest. In this range, 
the interface states can easily communicate with the minority carrier 
band but the minority carrier band has no connection with the bulk. 
Thus, the main effect minority carriers have in the equivalent circuit 
is to tie all R n ,,C, together as shown in Fig. 9(d), thereby eliminating 
the time constant dispersion. 

In the region around mid-gap where majority and minority carrier 
densities at the silicon surface are comparable, no simplification is 
possible and the equivalent circuit of Fig. 9(c) holds. This region is 
only a few kT/q wide and no analysis was attempted of measurements 
in this region. 

V. MEASUREMENT APPARATUS, REPRODUCIBILITY AND SAMPLE DRIFT 

5.1 Measurement Apparatus 

The measurement of the equivalent parallel conductance of an MIS 
capacitor for the purpose of extracting interface state properties re- 
quires consideration of the following points which are much more im- 
portant in conductance measurements than in capacitance measure- 
ments alone. 

(i) Harmonics of the signal frequency due to the non-linearity of the 
charge-voltage characteristic of the MIS capacitor can give rise to a 
conductance. To insure that this is negligible and the conductance 
measured is due to interface state transitions only, the maximum 
swing of surface potential caused by the applied ac signal should be 
less than kT/q volts. 

(ii) The conductance measurements reported here are in the range 
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between a nmho and about 10 /i mhos which corresponds to very small 
loss angles. Such small conductance values require that current leakage 
paths along the air-oxide interface be minimized. 

(in) To verify the theory developed in the previous section, capaci- 
tance, equivalent parallel conductance, and bias because of the ex- 
ponential dependence in (22) must be measured very accurately. 

The slices with the completed MIS capacitors were placed on a brass 
pedestal capped with platinum in a light-tight grounded box for 
measurement. A steady stream of dry nitrogen was supplied to the box. 
This keeps the oxide surface dry thereby preventing the formation of a 
current path for the ac signal along the oxide surface from the field 
plate to the back contact. It is especially important when making con- 
ductance measurements that such current paths be eliminated. 

To avoid scratching the field plate and thereby reducing its area, a 
gold wire, 125 /x in diameter, mounted on a micromanipulator was 
used to make contact to the field plate. 

For high sensitivity, it is necessary to measure accurately con- 
ductance when the loss angle is very small. The most suitable instru- 
ment for doing this over a wide range of frequencies is a capacitance 
bridge. The General Radio 1615-A capacitance bridge was found 
satisfactory over the frequency range from 50 Hz to 20 kHz and the 
Boonton 75-C capacitance bridge from 5 kHz to 500 kHz. Three term- 
inal capacitance measurements were made with these bridges. Fig. 10 
is a block diagram showing the arrangement. 

There are three parameters which require special attention in these 
measurements. They are (i) frequency, (u) signal amplitude, and (lit) 
dc bias. 

5.1.1 Frequency 

Frequency must be measured with a counter so that the precision 
in the calculated admittance of interface states and silicon is limited 
only by the precision of the bridges. The calibration marks on the 
oscillator dial in either the Boonton bridge or the hp signal generator 
are not accurate enough to obtain frequency ("see Section VI). 

5.1.2 Signal Amplitude 

To keep harmonics of the signal frequency from giving rise to a 
spurious conductance, only signals of small amplitude can be ap- 
plied. The small signal range is determined by experiment. Fig. 11 is a 
plot of measured capacitance and equivalent parallel conductance 
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Fig. 10 — Block diagram of measuring apparatus. 

normalized to their largest values as functions of peak-to-peak signal 
amplitude at a bias corresponding to peak conductance at 100 Hz, 
measured at a point where the slope of the MIS capacitance is large. 
There is a range of low signal amplitude at which both capacitance and 
equivalent parallel conductance are independent of signal amplitude. 
This is the experimentally determined small signal range. The maxi- 
mum amplitude of the ac silicon surface potential produced by this 
signal is about two fifths the applied signal amplitude. 

To see this, note that the applied signal divides between the oxide 
layer and the silicon. Expressing this mathematically, substitute (12) 
into (15): 



fan = 



., n , p>c d + y. ' 



(45) 



For a fixed signal amplitude, the amplitude of the silicon surface po- 
tential varies with bias because Cd and Y„ are functions of bias. 
Maximum silicon surface potential amplitude occurs when the Fermi 
level is near mid-gap because that is where Cd and | Y a | are smallest. 
The reason j Y s I is small is that the majority carrier density is small 
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making the interface state time constant long. Thus, to a good ap- 
proximation I o)(C ox +C D ) I > I Y 8 I so that (45) becomes 



8f. 

8v n 



(' 



C BI + C D 



(46) 



Typical values are : C ox = 5.73 X 10~ 8 farads/cm 2 and C D = 8.60 X 10 _s 
farads/cm 2 which gives 8\{/,/8v a = 2/5. 

At the maximum sensitivity setting of the null amplifier, null is 
found in the small signal range at a null meter reading of £t;l/2 j*A on 
the Boonton bridge and £^3/4 /*V on the General Radio bridge. These 
are the null readings determined essentially by the noise in the sys- 
tem. Therefore, harmonics of the signal frequency arising from the 
nonlinearity of the charge-voltage characteristics of the MIS capacitor 
are unimportant. However, as signal amplitude is increased beyond 
this range, nulls are obtained at successively higher null meter read- 
ings. This means that the amplitude of the harmonics has become 
large enough to get through the tuned null amplifier. In this range of 
signal amplitude, the conductance peak decreases in magnitude and the 
conductance curve is broadened. 

Operating in the small signal range entails a reduction in bridge 
sensitivity. Therefore, optimum signal amplitude is at the high end of 
the small signal range. The small signal range is determined by the 
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Fig. 11 — (a) Normalized capacitance and (b) normalized equivalent parallel 
conductance at fixed bias. The maximum value of capacitance is 367.8 pF and 
conductance 58.8 n mhos. Frequency is 100 Hz. Bias corresponds to peak con- 
ductance at this frequency. 
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oxide thickness and silicon resistivity. This range extends to higher 
signal amplitudes as the oxide is made thicker or silicon resistivity 
reduced. For oxide thicknesses between 500 A and 700 A and silicon 
resistivities around 1 ohm cm which is typical of the samples meas- 
ured in this experiment, optimum peak-to-peak signal amplitude is 
found to be 50 mV. For other oxide thicknesses and silicon resistivities, 
optimum peak-to-peak signal amplitude can be found by repeating the 
measurements shown in Fig. 11. An external attenuator is needed on 
the hp signal generator to get an amplitude of 50 mV. Signal amplitude 
is measured with an oscillosocope and adjusted whenever frequency is 
changed. The maximum peak to peak silicon surface potential with this 
signal is 2/5 X 50 = 20 mV which is less than kT/q at 300°K and 
therefore small signal. 

5.1.3 DC Bias 

Because interface state time constants are an exponential function 
of silicon surface potential, dc bias must be measured to at least three 
decimal places. The Cimron 6900A digital voltmeter is more than 
adequate for this purpose. 

DC bias is supplied to the MIS capacitor through a voltage divider 
in each bridge. Therefore, it must be measured across the sample 
rather than at the supply terminals. The voltmeter puts an additional 
capacitance and conductance across the sample which are read by the 
bridge. To avoid the correction this entails, the voltmeter is switched 
out of the circuit before the bridge is balanced. Because the input 
resistance of the voltmeter no longer shunts the MIS capacitor, the 
current through the voltage divider in the bridge decreases. Therefore, 
the voltage across the MIS capacitor is slightly larger than that meas- 
ured. For the Cimron 6900A digital voltmeter which has an input 
impedance when balanced of 10° ohms, this is not a problem. However, 
for a digital voltmeter having an input impedance of a few megohms, 
this correction becomes apparent whenever the resistance values of the 
internal bridge voltage divider are changed. These changes occur when 
switching the conductance multiplier on the Boonton bridge, the 
conductance to loss factors switch on the General Radio bridge or in 
switching from one bridge to another. To eliminate this error, a Gen- 
eral Radio 1230A electrometer which has an input resistance of 10 10 
ohms with the resistance selector set at oo can be used to measure 
bias voltage across the sample. The output of the electrometer is then 
read by the digital voltmeter. 

To permit control of the bias voltage to within a fraction of a 
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millivolt, a voltage divider consisting of a bank of fixed resistors in 
series with a 10-turn helipot is used in conjunction with a heavy duty 
45-volt battery in a shielded box. 

5.2 Reproducibility of Data and Sample Drift 

Sample preparation is designed to eliminate drift under negative 
bias at room temperature only. The admittance-voltage characteristics 
in the range of interest for both n- and p-type samples are located in 
the negative bias range because of the presence of a fixed residual 
positive charge in the oxide. Therefore, n-type samples are in strong 
accumulation and p-type in strong inversion at positive bias. No con- 
ductance measurements are made at positive bias because conductance 
due to interface states is too low to detect in strong accumulation or 
inversion. 

Reproducibility of the measurements on a given MIS capacitor was 
assured as follows. Measurements are repeated after room temperature 
bias aging at high negative bias, after the lapse of several days, and by 
hysteresis measurements where the admittance-voltage characteristics 
are measured first with increasing bias and then with decreasing bias. 
Only those samples where the measurements were found to be repro- 
ducible within experimental error are used. Interface state properties 
for such capacitors on a few different slices are found to be similar. 

VI. TECHNIQUE FOR ANALYZING THE DATA 

6.1 Properties of Equivalent Parallel Conductance 

The variation of G p /<a with bias or frequency can be most easily seen 
from (17) and (22) for a single time constant. In (39) derived from 
the statistical model which fits the measurements exactly in the 
depletion-region, (see Fig. 21) the variation of G p /u> with bias and 
frequency is obscured by the mathematical complexity. Equation (17) 
is adequate for illustrative purposes because of its simplicity and the 
fact that its variation with bias and frequency is qualitatively the 
same as (39) in the depletion region. Equation (17) fits the data 
exactly in the weak inversion region (see Fig. 25). G p or G p /<a can be 
obtained from measurements as a function of bias with frequency 
as parameter or as a function of frequency with bias as parame- 
ter. G p or G p /oi will go through a peak as a function of bias at 
each frequency as seen from (17) and (22). For a continuum of states, 
this peak occurs when the response time of the interface states which is 
being varied by the bias equals the period of the applied signal 
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frequency multiplied by a constant. When the frequency is increased, 
this peak shifts to a bias value corresponding to a surface potential 
nearer to flat bands. Fig. 1 shows capacitance and equivalent parallel 
conductance measured at 5 kHz and 100 kHz in a wet oxide. This 
figure illustrates how measured equivalent parallel conductance varies 
with bias and frequency. Fig. 1 also shows how inaccurate the capaci- 
tance method for extracting information about interface states is 
compared to the conductance method. It can be noted from Fig. 1 that 

(i) Maximum capacitance dispersion in this frequency range is 
about 14 percent while equivalent parallel conductance changes 
more than an order of magnitude. 

(ii) The conductance peak moves towards a bias voltage closer to 
flat band bias which in Fig. 1 is -2.6 volts as frequency is increased. 
Between 50 Hz and 500 kHz, conductance does not peak at flat band 
bias in the samples measured.* The peak shifts over the bias range 
from weak inversion to depletion in this frequency range. 

(Hi) Capacitance and equivalent parallel conductance are related to 
each other by the Kronig-Kramers relations. Rather than going 
through a mathematical proof of this, we note that (16), (17), (39), 
and (40) have the same form as the Debye equations. 10 ' 20 Because tke 
Debye equations satisfy the Kronig-Kramers relations, 27 we infer that 
(16), (17), (39), and (40) do also. 

Gp/o> as a function of frequency will go through a peak at each value 
of bias but in this case. G p will not. It can be seen from (17) that only 
G p /<a is a symmetric in wr when frequency is the variable while both 
G p and G p /o> are symmetric in or when bias is the variable. 

Gp/w vs bias (with frequency as parameter) will be spread over a 
bias range determined by the time constant dispersion and the density 
of interface states. The density of states spreads the G p /a> vs bias curve 
in the same way that it spreads the high-frequency capacitance vs 
bias curve. However, G p /<o vs w (with bias as parameter) will be spread 
over a frequency range determined only by the time constanct disper- 
sion. Because time constant dispersion determines the interface state 
admittance, G p /<o vs w curves are the more useful of the two for the 
purposes of this work. 

6.2 Extraction of G P /o> 

It is necessary first to get G p /a> and C p from the measurements. C p 
and G p are shown in Fig. 8(b) to be the capacitance and equivalent 

* To get the conductance to peak at flat band bias, frequencies higher than 500 
kHz have to be applied. 
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parallel conductance of the portion of the equivalent circuit consisting 
of C D in parallel with the distributed network representing the inter- 
face states. The admittance measured across the terminals x-x in 
Fig. 8(a) is G m + jwC,„. Converting this to an impedance, subtracting 
away the reactance of C ox , and converting back to an admittance, we 
get 



Gj, o>ClG m (Gl + o?Cl) 

<o ■ JClGl + [o 2 C m (C 0Z - C m ) - Glf 



(47)' 



and: 



r C OI (G 2 m + a?Cj)WC m {C oz - C.) - Gl] 

° p aC 2 oz G 2 m + [co 2 C m (C OI - C m ) - Gl] 2 ' (48; 

The measured G p /u> from (47) is related to N„ and r„, by (39) and C p 
from (48) is related to C D , N Bg and r m by (40) in the depletion region. 

It should be noted that the circuit across the terminals x-x in 
Fig. 8(a) has a shorter time constant than that of just the interface 
state branch of the circuit because of C ox . Therefore, G p /a from (47) 
will peak at a lower frequency than G p /o> measured across terminals 
x-x when equivalent parallel conductance is measured as a function 
of frequency with bias as parameter. If equivalent parallel conductance 
is measured as a function of bias with frequency as parameter, G p /a> 
from (47) will peak at a bias close to flat bands for the same reason 
as can be seen from (22). 

To minimize errors introduced by correcting for C ox , the oxide 
should be made as thin as practicable and frequency measured ac- 
curately with a counter as described in Section V. 

6.3 Bias vs Surface Potential 

From (22) , there is an exponential relation between interface state 
time constant and surface potential. Therefore, bias is measured with 
a digital voltmeter as described in Section V and the relation between 
bias and surface potential is determined as accurately as possible. It 
is necessary to do this to prove that (22) fits the measurements and 
to extract c„ and c p accurately. One of the most accurate methods for 
obtaining the relation between bias and surface potential has been 
developed by C. N. Berglund. 28 In this method, capacitance is meas- 
ured as a function of bias at such a low frequency that the interface 

* These enuations were evaluated using the measured values of u, C„ x , C m , and 
Gm on an IBM 1620 computer. 
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states are in equilibrium with the silicon. This means that the loss 
associated with the capture and emission of carriers by the interface 
states is negligible. 

Bias vs surface potential is obtained in two steps from the measured 
data. First, surface potential is found to within an additive constant 
at each bias. Then, the additive constant is found. 

6.3.1 Surface Potential within an Additive Constant 

Low-frequency capacitance C om measured at bias v is 

C om = dQ T /dv . (49) 

Differentiating (13) with respect to v 

1 = d^./dv + (1/C„) dQr/dv . (50) 

From (49) and (50) 

d*,/dv = 1 - C om /C oz . (51) 

The surface potential at any bias v, A is found by integrating (51) 

yp.ivox) = f " (1 - C B JC BX ) dv + A, (52) 

■* Oca 

Where A is an additive constant and v o2 is a bias in strong accumula- 
tion. 

Thus, surface potential at each bias can be obtained to within an 
additive constant directly from measured data by an integration. The 
fact that an integration is carried out using directly measured quanti- 
ties makes this method accurate. 

It is found that in samples oxidized in steam as described in Section 
II, 50 Hz is a low enough frequency to measure C om vs v and use (52) 
with negligible error. This is seen by the observation that the loss 
tangent is negligibly small in the bias range investigated here below 
about 80 Hz so that dispersion in the capacitance-voltage characteris- 
tic virtually disappears in this frequency range. 

Another important consideration is that (52) will not give an 
accurate result if there are any gross nonuniformities. Gross nonuni- 
formities occur when there are macroscopic areas under the field plate 
in which the charge density is significantly different from adjacent, 
areas. Whether or not there are such nonuniformities present in a 
sample can be determined by evaluating (52) using C om measured 
from strong accumulation to strong inversion. The resulting value of 
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fa will be slightly less than the bandgap,* if there are no gross nonuni- 
formities and considerably higher if there are. 

The random distribution of built-in charges and charged interface 
states does not affect the accuracy of (52). Because of the statistical 
fluctuations of surface potential caused by this random charge distri- 
bution, bias is really found as a function of mean surface potential. 

6.3.2 Determination of Additive Constant 

The additive constant A in (52) and the doping density is calculated 
next. For this purpose, C D must be found as a function of bias in the 
depletion region. Then CJ,* is plotted as a function of \J/„ — kT/q where 



>.. = f " (1 - C JC oz ) dv a . 



GJu from (47) peaks at bias values within the depletion region at 
frequencies above a few hundred Hz. C D is found from admittance 
measured at these frequencies using (48) to get first C p . Then from (40) 

C p = C D + C,(u>, u.), (53) 

where C, (co, w.) is the equivalent capacitance of the distributed network 
representing the interface states in Fig. 8(a) at the frequency and bias 
corresponding to (G v /u>) mn of the measured curve obtained using (47). 
It is shown in Appendix C that C,(w, u.) = 2(G p /co) max where (G r p /co) max 
is the measured peak. Then, using (53) we get 

C D = C p - 2((? p / w ) max . (54) 

In the depletion region except near flat bands, C D can be written 
in the following form using the approximate relation between C D 
and u, 16 

C~ D 2 = 2(*7, - l)/qe ti N A (3, (55) 

where e,, is the dielectric permittivity of silicon. Thus, a plot of C^ 2 
measured far from flat bands in the depletion region vs $. t — kT/q 
measured in the same region will yield a straight line as shown in Fig. 12. 
The intercept of this line with the abcissa when it is extrapolated to 
Ct> 2 = gives A. From this, bias is found as a function of surface po- 
tential. Finally, the slope of the line in Fig. 12 gives the doping density 
from (55). 



*This is particularly (ho rase for samples having high doping density such as 
those used in this work. 
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6.4 N„ and r m in Depletion 

Equation (39) for G p /u from the statistical model is found to fit 
the measured GJu vs frequency curves accurately at each bias (see 
Section VII). To prove this, (39) is numerically integrated on an IBM 
7094 computer using the trapezoidal rule.* u, corresponding to the 
particular bias at which GJu vs frequency is measured is obtained from 
the relation between bias and surface potential just found. The char- 
acteristic area a is the only parameter which is varied to obtain the 
best fit to the measured curve. JV„ can be obtained independently. 
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Fig. 12 — 1/Cd 2 vs \[>. e —kT/q where C D in farads is obtained from measure- 
ments in the depletion region of an n-type sample (N D = 1.15 X 10 1B cm -3 ) and 
t.c from the integral in (52). Field plate diameter is 3.8 XlO -2 cm. 

The parameter a is not arbitrary. It is found that a* has about twice 
the value and the same bias dependence as the silicon depletion layer 
width, (see Fig. 22). 

Fig. 13(a) shows G p /u found by integrating (39) plotted as a function 
of ojT m with N„ arbitrarily chosen as 1 X 10 11 cm -2 eV" 1 and a at a 
bias midway between flat bands and mid-gap. For comparison, Fig. 
13(b) is a plot of (17) and Fig. 13(c) is a plot of (20) vs uT m for the 
same value of N„ . From Fig. 13(a), it is seen that the maximum value 
of GJw from (39) occurs when wr m = 2.5. This was found to hold 

* To do this, limits of integration in (39) must be finite. Because the fluctuations 
of surface charge are small, integrating from u. = — 2u B to u. = 2Ua will give 
a good approximation to the integral from — oo to oo. 
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Fig. 13 — Calculated G P /w vs log u>T m . (a) Plot of (39) incorporating time 
constant broadening from statistical fluctuations of surface potential as well as 
continuum of states, a == 2.5 X 10" 10 cm 2 and impurity density = 2.1 X W cm" 3 are 
used, (b) Plot of (20) calculated by Lehovec in Ref. (19) for a continuum of 
constant density and capture cross section. The same interface state density 
arbitrarily chosen as 1 X 10 u cm^-eV" 1 is used in calculating all the curves. 



at all values of bias measured. The condition cor m = 2.5 is used to 
find T m from the frequency at which the measured G v /a) vs frequency 
curve goes through a maximum at each bias. The r m found this way 
corresponds to the mean surface potential u, at each bias. The condition 
o)T m = 2.5 can be used to get r m when G p /u is measured as a function 
of bias with frequency as parameter. In this case, u is given by the 
fixed frequency but the r m obtained now corresponds to the bias at 
peak G p /(i). 

Interface state density vs bias is found by numerically integrating 
(39) with a)T m = 2.5, equating the result to the measured peak G p /w 
value, and solving for iV,, at each bias. The condition wr m = 2.5 is 
equivalent to setting u„ = u„ — u s + In 2.5 in (39). Doing this, (39) 
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becomes 



and 



(^) = iqN..[2wtf + ■i)]"**— (56) 

\W /max 

/ max = f UB exp [-(a + !/„)] In [1 + exp (2y m )] du, , (57) 

J-2ufl 



where 

y m = u B — u, + In 2 .5 and z = (w 8 — u,) 2 /2(<r 2 , + <r|) . 

Equation (56) and (57) can be evaluated once a is known. The 
simplest way to find an approximate value for a is to use the experi- 
mentally established relation between a* and W where W is the deple- 
tion layer width. W can be obtained from the relation 

W = ^-, (58) 

where C D is obtained from (54). 
Solving (56) for N„, 

N, s = 2(^) [2t(<t* + iMty.-rV 1 . (59) 

\« /max 

Substituting (57) into (59) and using values of a from Fig. 22, (59) 
was evaluated to find N„ as a function of the position of the Fermi 
level in the bandgap of the silicon. 

6.5 N M , and r m in Weak Inversion 

This is a particularly simple case because G p /u vs frequency is char- 
acterized by a single time constant at each bias. G p /o) vs frequency 
or bias curves are obtained using (47) as before. Fig. 6 and (17) apply 
so that the condition for maximum G p /a is &>T m = 1 and the value 
of (G p /u) max is C./2. 

6.6 Capture Probabilities and Cross Sections 

Log r m is plotted as a function of u B — u, from measurements in 
the depletion region to see if the data fit (22). It is found this way 
that the measurements fit (22) quite well for both holes and electrons 
(see Section VII). The fact that the experimental data fit (22) means 
that c n and c p are independent of bandgap energy. Therefore, c n and 
c p are not affected by statistical fluctuations of surface potential and 
are the true capture probabilities of the interface states. Capture 
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cross section is related to capture probability by the relation c„ = v<r H 
for electrons where v is the average thermal velocity of electrons. 
c„ and c p are calculated from experimental plots of (22) and <t„ and a v 
found using v = 10 7 cm/sec for both holes and electrons. 



VII. EXPERIMENTAL RESULTS 



7.1 Depletion Region 



7.1.1 Sleam-Grown Oxides on [111] Orientation 

Fig. 14 shows log Tm vs u B — u, calculated from the measurements 
as explained in Section VI. The values plotted correspond to peaks 




Fig. 14 — Log T m vs ub — u. — showing fit of experimentally obtained points to 
(22) for both holes and electrons using steam grown oxide and [111] orientation. 
Points in lower half of the gap were measured on a p-tvpe sample. Field plate 
diameter is 1.2 x 10" 1 cm, acceptor density 2.08 X 10" cnr 3 , u B = 14.1 and C., = 
5.74 X 10~" farads/cm 1 . Points in upper half of the gap were measured on an n-tvpe 
sample. Field plate diameter 3.8 X 10 -2 cm, donor density 1.15 x 10 18 cm" 8 u B = 
13.5, and C„, = 5.08 x lO^farads/cm 2 . 
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of the G p /co curves obtained from measurements made at 300°K. The 
condition ur m = 2.5 is used to get r m . Because only majority carrier 
transitions give rise to measured conductance, curves in the lower half 
of the silicon bandgap were obtained from measurements on p-type 
and in the upper half on n-type. It is seen that these curves fit (22) 
so that capture cross sections are obtained simply by extrapolating the 
curves to u B — u, = where (22) becomes r mo = (wi,<r) _1 . This is 
then solved for a to get capture cross section. Capture cross sections 
obtained this way from Fig. 14 are 



<t p (holes): 
cr„ (electrons): 
o n (electrons): 



1.7 X 10 _1B cm 2 
5.9 X 10~ 16 cm 2 



from curve (a) 
from curve (b) 
from curve (c). 



In all cases, v = 10 7 cm/sec and n { = 1.55 X 10 10 cm" 3 have been used. 

Fig. 15 shows the corresponding values of N„ as a function of the 
average position of the Fermi level E F with respect to mid-gap potential 
at the silicon surface. Measurements in the lower half of the gap were 
made on p-type and in the upper half on n-type. 

Fig. 16 illustrates the temperature dependence of r. Curve (a) is 
calculated from measurements made at 204 °K while curve (b) is just 
a replot of curve (a) from Fig. 14 for comparison. 
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Fig. 15 — N.. from (59) vs E F — E t measured in units of kT/q for same two 
samples in Fig. 14. Flat band is 14.1 in lower half of the gap and 13.5 in upper 
half. E r is position of the Fermi level and Ei is mid-gap. 
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Q05 0.10 0.15 0.20 0.25 0.30 0.35 

+S IN VOLTS 

Fig. 16 — Log T m vs #, for p-type sample of Fig. 14: (a) measured at 204 °K and 
(b) at 300 °K. (b) Replot of (a) in Fig. 14. Curves show that experimentally obtained 
points satisfy temperature dependence of (22). (a) is the only one in this experiment 
measured at a temperature other than 300 °K. 



7.1.2 Steam-Grown Oxides on [100] Orientation 

Fig. 17 shows log r m vs u B — u, for the [100] orientation calculated 
from measurements made at 300°K. Capture cross sections are 



g v (holes): 2.0 X 10 -10 cm 2 

a v (holes): 4.0 X 10" 16 cm 2 

a n (electrons): 1.2 X 10" 15 



cm 



from curve (a) 
from curve (b) 
from curve (c). 



Fig. 18 shows the corresponding values of N„ in a plot similar to 
Fig. 15. 
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Fig. 17 — Log r m vs u B — Us showing fit of experimentally obtained points to 
(22) for both holes and electrons using steam-grown oxide and [100] orientation. 
Points in lower half of gap were measured on a p-type sample. Field plate diame- 
ter is 1.37 X 10- 1 cm, acceptor density 1.19 X 10 18 cm" 8 , u B = 13.5, and C » = 4.94 X 
lO^farads/cm 3 . Points in upper half of gap were measured on an n-t.ype sample. 
Field plate diameter is 1.43 X 10" 1 cm, donor density 1.44 X 10" cm" 3 , u B = 13.7, 
and Cot = 4.24 X lO^farads/cnr. 



7.1.3 [111] Orientation, Partly-Dried Oxides 

Fig. 19 shows log r m vs u B — u. for oxides initially grown in steam 
and partially dried by heating in dry N 2 as described in Section III. 
Only n-type was analyzed in detail, so that we have 

<r„ (electrons): 4.6 X 10" 16 cm 2 from Fig. 19. 

Curve (a) in Fig. 20 shows the corresponding values of N„ . Curve (b) 
in Fig. 20 is a replot from Fig. 15 for comparison. 
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Fig. 18 — N,. from (59) vs Ep — E ( for same two samples in Fig. 17. Flat band 
is 13.5 in lower half of gap and 13.7 in upper half. 
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Fig. 19 — Log T m vs u B — u. after partial drying of n-type sample in Fig. 4. 
Measured points still fit (22). Donor density is 1.4 X KPcm -8 , u B = 13.7, and 
Co* = 5.96 X lO^farads/cm*. 
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Fig. 20 — N.. from (59) vs E F - E ( . (a) n-type sample in Fig. 19 and (b) replot 
of n-type sample from Fig. 15 for comparison. For (a), flat band is 13.7 and for 
(b), flat band is 13.5. 



7.1.4 Range of N„ 

As seen from Figs. 15 and 18, N„ in steam-grown oxides is in the 
10" cm"~ 2 -eV _1 range. Partly drying increases N„ as seen from Fig. 20. 
Complete drying or growing the oxide in dry oxygen at 1000°C results 
in JV„ values in the 10 12 cnT 2 -eV _1 range or about an order of mag- 
nitude higher than for steam-grown oxides. 

N„ can be reduced below the values for steam-grown oxides by an 
order of magnitude by a method described by Balk. 29 First, an aluminum 
film serveral thousand A thick is evaporated over the entire surface 
of a steam-grown oxide. The sample is then annealed in H 2 or N 2 at 
350°C for 1/2 hour. After annealing, field plates are produced using 
photoresist techniques and chemically etching off the unwanted alum- 
inum. This annealing process is found to reduce N„ from the 
10" cm" 2 -eV _1 range characteristic of steam-grown oxides to the 
10 10 cm" 2 -eV _1 range. Thus, N„ can be varied two orders of magnitude 
from the 10 12 cm~ 2 -eV -1 range characteristic of dry oxides to the 
10 10 cm~ 2 -eV _1 range. Capture cross section appears to remain ap- 
proximately the same throughout. 
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7.1.5 Fit to Statistical Model 

The circles in Fig. 21 represent points calculated from measurements 
on the [111] p-type sample with a steam-grown oxide in Figs. 14 and 15. 
These points were measured with frequency as variable and bias fixed 
at —1.9 volts which corresponds to u B — u. = 6.04. The best fit to 
these points calculated from (39) is the solid curve in Fig. 21. This 
fit was obtained by substituting (59) for N„ in (39) and using (G r P /co) max 
from the data. The value of a required for this fit was found to be 
2.5 X 10" l0 cm 2 . 

Fig. 22 shows a* vs space-charge width W. The values of a are ob- 
tained by finding the best fit of (39) to G p /u vs frequency curves taken 
at different bias settings on the same sample. Space-charge width is 
calculated from the depletion layer capacitance obtained as described 
in Section VI from the measured capacitance and equivalent parallel 
conductance. 



7.2 Weak Inversion 

The GJw curves obtained after correcting for C oz will peak in the 
weak inversion range at signal frequencies below about 1 kHz. Curves 
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Fig. 21 — Gp'/w ys frequency showing how experimental points fit (39). An a 
of 2.5 X 10" 10 cm 2 is used in (39) to get this fit. G P ' is equivalent parallel con- 
ductance in mhos derived from the measurements. Experimental points are 
measured at a bias of —1.9 volts 0?, = 8.1) in depletion region on the p-type 
sample of Fig. 14. 
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Fig. 22 — a J vs space-charge width where a is obtained by finding the best fit 
of (39) to experimental points taken on p-type sample in Fig. 14. a* can have 
any value within vertical bars through each point. Bars indicate precision with 
which (39) can be fitted to measured G P /« vs frequency points. No reason why 
the point at the largest W does not lie on the straight line is known. 



of capacitance and equivalent parallel conductance vs field plate bias 
measured at 50 Hz on a [111], n-type sample with steam grown oxide 
are shown in Fig. 23. An n-type sample was used to avoid the masking 
effects of lateral ac current flow. 10 The capacitance peak below mid-gap 
persists down to 1.6 Hz, the lowest frequency used, and disappears 
at frequencies above 500 Hz. Above 500 Hz, the capacitance curve 
becomes the usual high-frequency curve. 6 The low-frequency capaci- 
tance dispersion is shown in Fig. 24. 

GJu vs frequency curves in this range can be fitted with a single 
time constant as shown in Fig. 25. The circles represent points calculated 
from measurements on an n-type sample with steam-grown oxide on 
the [111] orientation. These points were measured with frequency as 
variable and bias fixed at —3.3 volts which corresponds to u B — 
u„ = —6.3 in weak inversion. The best fit to these points calculated 
from (17) is the solid curve in Fig. 25. The values: of t and C. used 
in (17) to obtain this fit are 1.7 X 10" 3 seconds and 5.14 X 10" 8 
farads/ cm 2 , respectively. 

Fig. 26 is a plot of log r m vs u R — w. • These values are calculated 
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from measurements in the weak inversion region using cor m = 1 to 
get r m . Fig. 27 is a plot of the corresponding values of C, as a function 
of the position of the Fermi level with respect to mid-gap potential 
at the silicon surface. Values of C, are obtained from the peak of the 
GJu curve using ur m = 1 and (17). Values of N„ using N„ = C,/q 
are given on the right-hand ordinate axis of Fig. 27. Although these 
values are not exact, they give the correct order of magnitude of N„ . 

Fig. 28 shows a low-temperature drying experiment similar to those 
in Figs. 3 and 4 except the curves in Fig. 28 were measured at 50 Hz. 
Curves (a) and (b) of Fig. 28 are plots of capacitance and equivalent 
parallel conductance vs bias for a steam-grown oxide on n-type silicon. 
Curves (c) and (d) are the capacitance and equivalent parallel conduc- 
tance vs bias after heating in dry nitrogen at 350°C for 17 hours. 

Fig. 29 is a plot of capacitance and equivalent parallel conductance 
vs bias measured at 50 Hz on a sample having N„ in the 10 10 cm~ 2 - 
eV -1 range. The sample is n-type oriented in the [100] direction and 
the low value of N fS is produced by the annealing process described 
previously. 
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Fig. 23 — Capacitance and equivalent parallel conductance measured at 50 Hz 
VS field plate bias. Sample is the n-type one in Fig. 14. curves show capacitance 
and equivalent parallel conductance in the weak inversion region. 
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VIII. DISCUSSION 

The discussion is divided into three parts. The first part deals with 
the depletion region and the second with the weak inversion region. 
The third part discusses the limitations of the conductance technique. 

8.1 Depletion Region 

8.1.1 Steam-Grown Oxides on [111] Orientation 

It is seen from Fig. 14 that points calculated from the measurements 
fit (22) quite well for both n and p type samples. Time constant de- 
creases by a factor of 1/e for every kT/q increase in u B — u, . The 
exponential dependence of time constant on surface potential shown 
in Fig. 14 means that time constant is inversely proportional to majority 
carrier density at the silicon surface as seen from (21) and (22). This 




-2 



-3 -4 -5 

FIELD PLATE BIAS IN VOLTS 



Fig 24 Capacitance vs field plate bias with frequency as parameter measured 

with apparatus described in Ref. 28. on n-type sample of Fig. 14. Curves show 
capacitance dispersion in weak inversion. 
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Fig. 25 — Gp'/w vs frequency showing how experimental points fit (17). C, = 
5.14 X lO^farads/cm" and t = 1.7 X 10^ seconds are used in (17) to get this fit. 
G p ' is equivalent parallel conductance in mhos derived from measurements. Ex- 
perimental points are measured at a bias of —3.3 volts (w. = 19.8) in the weak 
inversion region on the n-type sample of Fig. 14. 

is the experimental evidence that majority carrier transitions make the 
dominant contribution to measured equivalent parallel conductance in 
the depletion region. 

The exponential dependence of time constant on surface potential 
is important evidence in addition to the drying experiments of Figs. 3 
and 4, that the dominant loss is due to transitions to and from interface 
states. This follows from the fact that no other relaxation processes in 
the system have time constants which depend exponentially on surface 
potential. To see this, we shall consider generation-recombination 
though impurity levels in the silicon near mid-gap, majority carrier 
transitions to and from impurity levels in the silicon near the Fermi 
level and the transit time of majority carriers across the silicon space- 
charge region. Bulk generation-recombination is insignificant because 
the Fermi level does not cross impurity levels near mid-gap for band 
bending in depletion as seen in Fig. 5(b). Therefore, generation -re- 
combination is taken to be zero in Figs. 6 and 8. The Fermi level can 
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Fig. 26 — Log Tm vs u B - u. for n-type sample in Fig. 14 in weak inversion 
region. Curve shows that equivalent parallel conductance in weak inversion is 
dependent on majority carrier density at the silicon surface. 

still cross any impurity levels located near it [see Fig. 5(b)]. Majority 
carrier transitions to and from levels at this crossover point will con- 
tribute to the loss. The location of the crossover point from the inter- 
face will vary with bias. However, the crossover point is always at 
the same energy with respect to mid-gap so that majority carrier 
density there will be constant and independent of bias. Therefore, the 
time constant for this process will be independent of bias. As time 
constant is not observed to be bias independent in depletion in the 
samples measured, majority carrier transitions to and from levels in 
the silicon bulk must make a negligible contribution to the loss com- 
pared to majority carrier transitions to and from interface states. 

Loss due to the transit of majority carriers across the space-charge 
layer is negligible. 11 ' 12 It can be shown using the relations given in 
Refs. 11 and 12 that the value of space-charge resistance for majority 
carriers when the silicon surface is intrinsic is four orders of magnitude 
less than majority carrier capture resistance in the samples measured. 

The fact that curves (a), (b), and (c) in Fig. 14 can be fitted so well 
by straight lines implies that capture cross section in these regions is 
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to a good approximation independent of energy in the bandgap. This 
independence is an experimental fact for which no satisfactory explana- 
tion exists at present. Because capture cross sections are independent 
of energy, they will be unaffected by statistical fluctuations of surface 
potential. The values of capture cross section obtained from the linear 
portions of Fig. 14 are therefore truly characteristic of the interface 
states present. Curves (b) and (c) in Fig. 14 measured on n-type are 
characterized by two types of interface states having different capture 
cross sections with an abrupt transition between them. The reason for 
this is not understood. 

Fig. 14 shows that N„ is slowly varying with energy. It is reasonable 
to conclude from this, and the observation, that the time constant 
varies monotonically with surface potential over the same portion of 
the band gap and that there is a continuum of interface states rather 
than a discrete level. It is also seen from Fig. 15 that N„ is in the 
10" cm" 2 -eV _1 range. This means that the capturing centers are spaced 
too far apart in the plane of the interface for the wavefunction of an 




E F -E^ (units of kT/aJ 

Fig. 27 — C. and estimated JV., vs E F — E t for n-type sample in Fig. 14 in the 
weak inversion region. Estimated N.. = C./q. C. increases monotonically with 
Er — Et. 
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electron in one center to overlap a neighboring center. Transitions from 
one center to another are therefore, highly improbable so that it may 
be assumed that we do not have a band but just a continuum of levels 
closely spaced in energy. Transitions will only occur between the 
majority carrier band and levels within a few kT/q of the Fermi level. 
Cascading transitions from center to center which could occur in a 
band would probably have a different dependence of time constant on 
surface potential than shown in Fig. 14. Therefore, cascading transitions 
were not considered in the theory developed in Section IV. 




-5 -4 -3 -2 
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Fig. 28 — (a) capacitance and (b) equivalent parallel conductance measured 
vs field plate bias before drying; (c) capacitance and (d) equivalent parallel 
conductance vs field plate bias of same capacitor after heating in dry nitrogen 
for 17 hours at 350° C. This is the same sample as Fig. 4 only measured at 50 Hz. 
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Fig. 29 — Capacitance and equivalent parallel conductance measured at 50 Hz 
vs field plate bias. Field plate diameter is 3.8 X 10""cm, silicon resistivity 1 ohm- 
cm, n-type, [100] oriented, and C, z = 2.8 x lO^farads/cnr. Oxide is grown in 
steam then H 3 annealed as described in text to produce low N,, . 

The integration of (4) and (26) was performed for the case where 
N„ and capture cross section do not vary appreciably over an energy 
interval of several kT/q. Figs. 14 and 15 show the experimental justifica- 
tion for this. 

In Fig. 16, log r m is plotted as a function of $, rather than u B — u, 
to illustrate more clearly the temperature dependence. An alternate 
and equivalent form of (22) which follows from (21) for p-type is 



va„N , 



exp (q+./hT). 



(60) 



Curve (a) in Fig. 14 fits (60) quite well with a slope q/kT where T = 
204°K. Curve (b) fits with the slope for 300°K. Thus, the temperature 
dependence expected from (60) is observed. This agrees with Shockley- 
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Read theory 18 which is used in Section IV. Capture cross section can 
be obtained from (60) by extrapolating the curves in Fig. 16 to $. = 0. 
Doing this, we get <r v = 8.8 X 10" 17 cm 2 from curve (a) and <r v = 
2.1 X 10~ 10 cm 2 from curve (b). The capture cross section from curve 
(b) is of course in agreement with a v from curve (a) in Fig. 14. However, 
it is 2.4 times larger than a v from curve (a) of Fig. 16. This is probably 
due to the error in finding the value of A from (52). Normally, A can 
be found only within 20 to 30 mV of its correct value by the technique 
of Section V. Part of this error is due to neglecting the difference be- 
tween hole and electron effective mass. It is considered likely that in 
reality capture cross section is nearly independent of temperature in this 
range. The apparent variation is ascribed instead to the errors in finding 
the value of A. 

8.1.2 Steam-Grown Oxides on [100] Orientation 

The same general conclusions can be reached for the [100] samples 
from Figs. 17 and 18 as for the [111] samples just discussed but there 
are two differences. The first difference is that the time constant curves 
in Fig. 14 for the upper half of the gap in [111] p-type are similar to 
the curves for the lower half of the gap in [100] p-type in Fig. 17. Also, 
the curves in the lower half of the gap in Fig. 14 for [111] p-type are 
similar to the curves in the upper half of the gap in [100] n-type in 
Fig. 17. The second difference is that interface state density in Fig. 15 
for the lower half of the gap in [111] p-type is about twice as great 
as in Fig. 18 for the lower half of the gap in [100] p-type. However, 
interface state density in Fig. 15 for the upper half of the gap in [111] 
n-type is about the same as in Fig. 18 for the upper half of the gap 
in [100] n-type. The reason for these differences is not understood. 

Built-in charge density in the [111] samples varies from sample to 
sample from 5 X 10 n cm" 2 to 9 X 10" cm -2 in both n and p-type. 
This is about twice as great as in the [100] samples confirming the 
observations of Balk. 29 Also, the increase of N„ towards the band 
edges as shown in Figs. 15 and 18 is consistent with the results of Gray 
and Brown. 30 

8.1.3 [111] Orientation, Partly-Dried Oxides 

Fig. 19 shows that interface state time constants fit (22) after partial 
drying. Capture cross section for n-type has not changed appreciably 
although two types of states were not found. This may be due to the 
fact that too small a range of surface potential was covered in the 
measurements. Because capture cross section did not change, partial 
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drying does not appear to introduce new types of states — but only 
increases the number of states already present. 

Comparison of curve (a) to curve (b) in Fig. 20 shows that interface 
state density for the partly-dried oxide is greater than for the wet 
oxide. This can be seen directly from Fig. 3 and 4 for p and n-type 
by the increase in the peak of the measured conductance with drying. 

Fig. 3 for p-type shows by the shift of the curves to higher negative 
bias that a net positive charge has been induced in the oxide by partial 
drying. Fig. 4 for n-type shows by the shift of the curves towards less 
negative bias that a net negative charge has been induced in the oxide 
by partial drying. No conclusions about these results will be drawn 
because it is not certain that other changes independent of interface 
state density have not occurred in these samples. 

The increase in interface state density with drying can explain the 
results observed in n + p junctions in a previous paper. 31 In Ref. 31, 
it was found that drying an oxide heavily contaminated with alkali 
metal ions drastically reduced the current flowing in an n channel 
along the surface of the p-type material. The increased interface state 
density caused by drying results in more channel electrons being 
captured by interface states. This markedly reduces field-effect mobility 
which in turn results in an increase in the sheet resistance of the inver- 
sion layer. Although the inversion layer may be terminated by a defect 
at the edge of the wafer which shorts it to the p substrate, the current 
flow will be reduced to a much lower level than before. 

8.1.4 Range of N., 

The range of N.. from 10 12 cm _2 -eV _1 for very dry oxides through 
10 11 cm~ 2 -eV _1 for wet oxides to 10 10 cm _2 -eV -1 for aluminum clad 
oxides annealed in H 2 or N 2 are empirical observations. It is not de- 
finitely known what these centers are. Therefore, an explanation can 
only be speculative. However, these results are qualitatively similar 
to previous work on "bare" surfaces which had at most a few mono- 
layers of oxide. 32 Exposure to water vapor always reduced N,. by a 
few orders of magnitude. 32 For the thick oxides studied in this work, 
water was most likely present in the form of OH groups which have a 
permanent dipole moment. An interesting model to explain the annihila- 
tion of fast states on "bare" surfaces has been proposed by Rzhanov. 33 
Assuming that each interface state also has a permanent dipole moment, 
the lowest energy configuration will occur at the interface when the 
dipole moment of an H 2 or OH group is aligned antiparallel with the 
dipole moment of the nearest interface state. This could reduce the 
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capture probability of the interface state 34 to such an extent that there 
is no measurable loss due to transitions involving it. The only loss 
measured would be due to transitions to and from the remaining un- 
paired states. For thick oxides, this fits the observed decrease in N„ 
with increasing wetness and the observation that measured capture 
cross section remains unchanged. This model is, however, only one of 
several capable of explaining the results. 

8.1.5 Statistical Model 

The time constant dispersion actually observed is much broader than 
expected for a continuum. This can be seen by comparing Fig. 21 to 
the curves in Fig. 13. Three mechanisms which could cause this broaden- 
ing are considered. They are: 

(i) Tunnelling of majority carriers to states randomly distributed 
into the oxide. 2 

(ii) Fluctuations of capture cross section over the plane of the 
interface. 

(Hi) Random fluctuations of surface potential over the plane of the 
interface. 

The tunnelling mechanism is unlikely to be the dominant cause of 
the time constant dispersion for the following reasons. The range of 
frequencies used in the measurements is from 50 Hz to 500 kHz. Thus, 
only those interface states having time constants shorter than several 
milliseconds will be measured. The time constant distribution for 
tunnelling would be 2 

r, - (n,^)" exp (2KJ-), (61) 

where <r„ is the true electron capture cross section of the states, 1/2K 
is the tunnelling decay constant, and £ is the distance into the oxide 
measured from the interface. For any reasonable value of 1/2K , such 
as 1 A given in Ref. 2, and reasonable values of <r , such as 10" 1B cm 2 
to 10" 16 cm 2 , it can be seen from (61) that states located more than 
a few k into the oxide would have time constants much too long to 
be measured in the frequency range between 50 Hz and 500 kHz. 
Hysteresis in the capacitance vs bias characteristic which can be ob- 
served by slowly varying bias from a low to a high value and then 
back is negligibly small. Hysteresis would be observed when there is 
an appreciable density of states deep in the oxide having time constants 
longer than the time it takes to vary the bias because the charging 
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and discharging of these states would lag behind the bias variation. 2 
The fact that no hysteresis is observed means that the density of states 
in the oxide having time constants of the order of minutes is negligible 
in the samples measured. Finally, it is shown in Appendix D that for 
any reasonable value of K , the shape of the GJu> vs log w curve cal- 
culated using (20) for the continuum of states and (61) for the tunnelling 
time constant distribution differs markedly from the observed G p /a> vs 
log co curve. 

Although capture cross section is observed to be independent of 
energy over large portions of the band gap from Fig. 14, capture cross 
section could depend upon the location of the state in the plane of 
the interface. The assumption that capture cross section fluctuates over 
this plane does not give a broad enough time constant dispersion to 
fit the observed G p /co vs frequency curves. 

Random fluctuations of surface potential over the plane of the inter- 
face is the model which fits best. One can see from (22) that small 
fluctuations of v, will cause large fluctuations in t„ . This broadens 
the time constant dispersion sufficiently to fit the experimental curves. 

Random fluctuations of surface potential can be caused by: 

(i) A random distribution of built-in charges and charged interface 
states over the plane of the interface. 

(it) A random distribution of ionized impurities in the silicon. 
(iit) Random fluctuations of oxide thickness over the plane of the 
interface.* 

Of these, the first two are the basis for deriving (39) in Section IV. 
The last was excluded because no dependence of time constant disper- 
sion on oxide thickness was observed between 500 and 2000 A in the 
samples measured. Fig. 21 shows the fit of (39) to experimental values 
of Gp/u measured in the frequency range between 50 Hz and 500 kHz 
at a mean surface potential of 8.1. Fig. 21 shows that (39) accurately 
fits the experimental points. The fit of (39) to the experimental G p /oj vs 
frequency points is equally accurate for all the other values of mean 
surface potential u, measured in depletion. The only independent 
parameter in (39) varied to obtain these fits is the characteristic area a. 
By fitting (39) to the experimental G p /u vs frequency points at each 
value of u. in this way, the values of a obtained are found to be related 
to space charge width W. This relation is shown in Fig. 22 where the 
length of the side of each square characteristic area a* is plotted against 



* This effect was pointed out to us by C. R. Crowell. 35 
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W. Fig. 22 shows that a* increases linearly with W and has a larger 
magnitude. To understand this relationship, consider a positive built-in 
charge or charged interface state located at a point in the plane of 
the interface. This charge is sandwiched between the plane of the field 
plate and the plane at the edge of the space-charge region a distance 
W away. An image charge is induced in each of these planes. The 
image charge in each plane induces another in the opposite plane and 
so on. Thus, the charge at the interface induces an infinite number 
of image charges in the field plate and space-charge edge planes. Let 
us now add a number of charges in the neighborhood of this charge 
in the plane of the interface and calculate the charge distribution 
induced by all the charges and their images in the plane at the edge 
of the space-charge region. The largest spacing between charges which 
results in a constant induced charge distribution in the plane at the 
edge of the space-charge region is a*. Thus, charges within a charac- 
teristic area are indistinguishable as separate charges which makes 
surface potential over a characteristic area uniform. The characteristic 
area is a measure of the granularity of the surface-charge distribution 
seen at the edge of the space-charge region. A calculation based on 
the image charge model shows that a* has the values in Fig. 22 and 
a linear dependence on W. The value of a can be obtained within a 
factor of 2 from the image charge model so that a is not completely 
arbitrary. 

We shall show next that the random distribution of built-in charges 
and charged interface states is the dominant cause of surface potential 
fluctuations in the samples measured. Rewriting (36) from Section IV 
for the standard deviation of surface potential caused by the random 
distribution of built-in charges and charged interface states gives 

a '- [W{u.)C ot + U W 

Curve (b) in Fig. 30 is the variance o* calculated from (62) as a function 
of W. Q calculated from (31) and (32) with u. = u. is found to be 
independent of bias from —1.8 volts to -2.2 volts. This bias range 
corresponds to mean surface potential from 10.3 to 3.8. Assuming a 
constant value of interface state density of 3 X 10" cm~ 2 -eV _1 in the 
lower half of the gap from Fig. 15, the change in total interface state 
density is 3 X 10" X (10.3 - 3.9) X 0.026 = 5.1 X 10 10 cm" 2 . This 
is negligibly small compared to Q/q = 9.8 X 10" cm 2 . Curve (b) 
in Fig. 30 is calculated using Q - 1.57 X 10" 7 coul/cm 2 , C„ = 5.7 X 
10~ 8 farads/cm 2 and the values of a from Fig. 22. 
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Fig. 30 — (a) total variance of surface potential {a," + as 2 ) caused by random 
distributions of built-in charges and charged interface states and ionized accep- 
tors in the silicon space-charge region, a. and a B are calculated from (62) and 
(63), respectively, (b) variance of surface potential <r, 2 caused by random distribu- 
tion of built-in charges and charged interface states alone, a. is calculated from 
(62). (c) variance of surface potential o-g 2 caused by random distribution of ionized 
acceptors in silicon space-charge region. a B is calculated from (63). (d) variance 
of surface potential a t s caused bv fluctuations of oxide thickness. a x is calculated 
from (64). 



For the standard deviation of surface potential caused by a random 
distribution of ionized acceptors in the silicon space-charge region, 
we rewrite (37) derived in paragraph B.l of Appendix B 



o- B = 



ggjffxgjfofllE - exp (-ti.)l 
2[W(fi.)C oz + 6.J 



(63) 



Curve (c) in Fig. 30 is the variance o-j calculated from (64) as a func- 



1116 THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUGUST 1967 

tion of W. The values used are: C ox = 5.7 X 10~ 8 farads/cm 2 and 
N A = 2.1 X 10 18 cm" 3 . 

Finally, for the standard deviation of surface potential caused by 
oxide thickness fluctuation, we rewrite (101), derived in paragraph B.2 
of Appendix B, 

_ qN A 0W 2 (u.)[l - exp (-<)] dx , , 

Curve (d) in Fig. 30 is the variance a\ calculated from (64) as a 
function of W. The fluctuations of oxide thickness arise from the random 
nature of the oxidation process. It is reasonable to take dx as about 
one silicon lattice constant or 5A. The values used in (64) to get curve 
(d) in Fig. 30 are: x = 595 A, N A = 2.1 X 10 lfl cm -3 , and dx = 5 A. 

Curve (a) in Fig. 30 is <r 2 + al vs W obtained by fitting (39) to the 
experimental (?„/« vs frequency points. Fig. 30 clearly shows that the 
random distribution of built-in charges and charged interface states is 
the dominant cause of surface potential fluctuations in the samples 
measured. It can be easily seen that a* should vary linearly with W as 
follows. Over the range of surface potential measured, the experimental 
GJoi vs frequency curves broaden very slightly with increasing surface 
potential. This is shown quantitatively in curve (a) of Fig. 30. Over 
the same range of surface potential, curve (c) in Fig. 30 is virtually 
independent of W. The dominant term a] in curve (a) is therefore, 
nearly independent of W as shown in curve (b) of Fig. 30. For a. con- 
stant, an inspection of (62) will show that for Q constant and WC ox < €„• 
as is the case in the samples measured, a* must vary linearly with W. 

The density of built-in charges and charged interface states is in the 
high 10 11 cm -2 range in the samples measured. Therefore, these charges 
are about 100 A apart on the average which is too far for coulombic 
interaction between them to be important. This makes the assumption 
that they are randomly distributed a reasonable one. The mean number 
of built-in charges and charged interface states N in a characteristic 
area a can be calculated from the relation _iV = aQ/q (see Section IV) 
using values of a obtained from Fig. 22. N is found to vary from 140 
charges at W = 6.9 X 10" 6 cm to 246 charges at W = 1.07 X 10~ 6 cm. 
For these values of N, there is negligible error in using the Gaussian 
approximation to the Poisson distribution (28) in deriving (39) in 
Section IV. The mean number of ionized acceptors calculated from the 
relation M = W 3 N A varies from 26 charges at W = 1.07 X 10" 6 cm 
to 7 charges at W = 6.9 X 10" 6 cm. For all values of M except perhaps 
the last, there is negligible error in using the Gaussian approximation 
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to the Poisson distribution (81) in paragraph B.l of Appendix B in 
deriving (39). Any error at W = 6.9 X 10~ 5 cm is not important to 
the final results because the contribution of al to the total variance 
is so small. The assumption of small fluctuations made in deriving (39) 
also seems to be reasonable because of the rather good fit obtained 
to the experimental points. 

One can decide from (62), (63), and (64) which influence will be 
dominant for samples other than the ones used here. In spite of the 
fact that the samples investigated here were dominated by random 
fluctuations of built-in charges and charged interface states, it should 
be borne in mind that the other influences may become important in 
other cases. When built-in charge and interface state density is low 
and the oxide relatively thick, then bulk doping fluctuations will 
dominate. When the oxide is thin and the impurity density is high, 
random fluctuations of oxide thickness become important. The general 
conclusion is that there will always be random fluctuation of surface 
potential in simiconductor-insulator interfaces. 

No data has been taken at flat bands or in the accumulation region. 
However, it is expected that the equivalent circuit in Fig. 8(a) will 
also be valid in this region with all the distributed time constants 
shifted to shorter values. 

8.2 Weak Inversion Region 

In this part, it will be shown that the loss in weak inversion can 
arise from generation-recombination either through bulk states or 
interface states and how the conductance technique can be used to 
distinguish between them. It will be shown also that the observed 
effects in this region can be explained with a continuum of states rather 
than a single level near a band edge. 

When the silicon surface becomes inverted, the Fermi level will 
cross impurity levels near mid-gap as seen from Fig. 5(c). Generation 
recombination through bulk states as well as through interface states 
can contribute significantly to the measured loss. These two cases can 
be distinguished by the bias dependence of the measured equivalent 
parallel conductance. When generation-recombination through interface 
states dominates the loss, the equivalent parallel conductance goes 
through a peak as a function of bias in weak inversion because time 
constant varies inversely with majority carrier density at the surface 
(see paragraph 4.6). This case is illustrated in Fig. 23 where equivalent 
parallel conductance goes through a peak in weak inversion and drops 
to a very low value in strong inversion. 
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When generation-recombination through bulk states dominates the 
loss, equivalent parallel conductance does not go through a peak as 
a function of bias in weak inversion. The reason for this is that the 
loss due to this process is bias independent. This follows from the fact 
that hole and electron densities at the crossover point [see Fig. 5(c)], 
where the Fermi level crosses impurity levels located near mid-gap, 
are independent of bias. This case is illustrated in Fig. 29. The slight 
indication of a peak in the conductance curve means that generation 
recombination through interface states makes a slight contribution to 
the loss. However, the primary effect is that of a bias independent 
loss. The variation of measured conductance with bias is mainly due 
to variation of capacitance with bias. 

Fig. 23 is measured on a sample having N„ in the 10" cm" 2 -eV _1 
range while Fig. 29 is measured on a sample having N„ in the 10'° 
cm _2 -eV _1 range. Thus, samples can be made in which one or the other 
process dominates the loss. It should be noted that the conductance 
technique can be used to investigate generation-recombination through 
interface states or bulk states in the weak inversion region by choosing 
an appropriate sample. The investigation of generation-recombination 
through bulk states is beyond the scope of this work so that most 
of the samples investigated have interface state loss dominant in weak 
inversion. Therefore, bulk generation-recombination is made zero in 
the equivalent circuits of Fig. 9. 

In general, samples having N„ in the 10 u cm" 2 -eV _1 range or higher 
will have interface state loss dominant in the weak inversion range. 
Additional evidence for this in such samples is shown in Fig. 28. The 
only difference between Fig. 28 and Fig. 4 is the low signal frequency 
which is used in Fig. 28 to reveal the behavior of equivalent parallel 
conductance in the weak inversion region. The conclusion is again 
drawn from Fig. 28 that interface state loss must be dominant as from 
Figs. 3 and 4 discussed previously in Section III. 

The striking feature of the experimental G p /u vs frequency curves 
at each bias in weak inversion is that they can be accurately fitted by 
a single time constant using (17). An example of this is shown in Fig. 25. 
Three possible explanations of this will be considered. 

(i) A single level state located at an energy in the silicon band gap 
corresponding to the bias at the capacitance peak in Fig. 23. 
(ii) A single level near the valence band edge. 

(Hi) A continuum of interface states over the silicon band gap. 

The first explanation can be ruled out by an examination of Figs. 24, 
26 and 27. Fig. 24 shows that the capacitance peak shifts to greater 
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negative bias with decreasing signal frequency. The peak must always 
occur at the same bias for a single level. In Figs. 26 and 27, the range 
of u B — u. and E F — E { corresponds to a range of bias which includes 
the capacitance peak in Fig. 23. Both r m in Fig. 26 and C, in Fig. 27 
increase monotonically through the corresponding capacitance peak in 
Fig. 23. For a single level state, C, would go through a peak and r m 
first would be constant for E F on the valence band side of the single 
level and then it would decrease exponentially with u, for E F on the 
conduction band side. 

The range of E F — E { and u B — u. in Figs. 26 and 27 is limited by 
the capacitance dispersion shown in Fig. 24. This means that (52) 
cannot be used without significant error too far beyond the capacitance 
peak (at reasonable frequencies). 

To rule out the possibility of a single level state near the valence 
band edge, G p /u vs frequency curves were obtained on a p-type sample 
This sample had steam-grown oxide on a 1 fi-cm crystal oriented in the 
[100] direction. The [100] direction has a lower built-in charge density 
so that masking effects due to lateral ac current flow into the inverted 
layer beyond the field plate are minimized. A single time constant was 
again found to be characteristic of the G P /u vs frequency curves in the 
weak inversion region. The experiment was repeated on an n-type 
sample oriented in the [100] direction having a steam-grown oxide with 
the same result. The fact that the same result was obtained in weak in- 
version on both n and p-type samples and the lack of evidence for a 
single level state in the same energy range in accumulation in these 
samples makes such a state as an explanation of the effect very unlikely. 
The observations can be best explained in terms of a continuum of 
states as done in paragraph 4.6. 

Evidence that this effect is dominated by the majority carrier density 
at the silicon interface is given by Fig. 26. Fig. 27 shows that C, in- 
creases so rapidly with E F that the approximation of a constant density 
of states over an interval of kT/q cannot be made. Thus, (18) cannot 
be simply integrated. The fact that C. is so rapidly varying is the 
reason why the curve in Fig. 26 does not have the same magnitude 
of slope as the curves in Fig. 14. This rapid variation of C, also explains 
the capacitance peak in Fig. 23. It can be seen that (16) will go through 
a peak if C, has a strong dependence on bias. 

Comparison of Fig. 27 with Fig. 15 shows that interface state density 
in the lower half of the silicon bandgap in n-type has a different magni- 
tude and dependence on position of Fermi level than in the lower half 
of the gap in p-type. The reason for this is not understood. 
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8.3 Limitations 0/ the Conductance Technique 

Interface state densities in the 10 9 cm~ 2 -eV _1 range have been meas- 
ured by this technique using the apparatus described in Section V. 
For a field plate diameter of 5 X 10" 2 cm, this corresponds to about 
2 X 10 8 states/eV. For a signal of 20 mV peak-to-peak, transitions 
involving only about 4 X 10 4 states can be detected. 

It is important to use as thin an oxide layer as practicable to get 
large values of capacitance and equivalent parallel conductance for a 
given N„ and to minimize the error in extracting C OI . 

The range of bandgap energy over which measurements were made 
in this work was limited by the frequency range 50 Hz to 500 kHz of 
the two bridges used. 

IX. SUMMARY ANI> CONCLUSIONS 

(1) Majority carrier transitions to and from interface states make 
the dominant contribution to the measured equivalent parallel con- 
ductance in the accumulation, depletion, and inversion regions. This 
fact makes it possible to use simplified equivalent circuits and extract 
interface state properties for each type of carrier independntly. 

(w) A continuum of Shockley-Read type interface state levels closely 
spaced in energy across the silicon band gap seems to be characteristic 
of the Si-SiOa interface. 

(m) Interface state densities in the 10" cm^-eV" 1 range are usual 
for "as-grown" steam oxides. 

(iv) Interface state density can be increased to the 10 1Z cm~ 2 -eV~' 
range by drying the oxide. 

(y) Interface state density can be reduced to the 10 10 cm~ 2 -eV -1 range 
by H 3 or N 2 annealing at low temperatures with an aluminum film 
deposited on a steam-grown oxide. 

\vi) Capture cross sections for electrons and holes are independent 
of energy over large portions of the silicon bandgap. They are also 
independent of temperature, the processes of drying or H 2 or N 3 an- 
nealing, and, except for subtle differences, crystal orientation. 

(vii) In depletion, a time constant dispersion is observed which is 
much broader than expected for a continuum. The broadened time 
constant dispersion is determined mainly by statistical fluctuations 
of surface potential. In the samples studied, a random distribution 
of built-in charges and charged interface states was found to be the 
primary cause of the surface potential fluctuations. This model provides 
an accurate fit to the experimental points and yields a characteristic 
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area whose linear dimension has about twice the magnitude and has 
the same bias dependence as the silicon space-charge width. Charges 
within a characteristic area cannot be distinguished as separate charges. 
Random fluctuations of surface potential can also arise from a random 
distribution of ionized impurities in the bulk silicon and by random 
fluctuations of oxide thickness. In samples where the built-in charge 
and charged interface state density is low or the oxide layer thin, 
these can replace the random distribution of built-in charges and 
charged interface states as the dominant causes of surface potential 
fluctuations. 

(viii) In weak inversion, the fast response time of minority carriers 
and the conductance of the inversion layer results in the observation 
of a single time constant at each bias. 

(ix) The interface state branch of the equivalent circuit in the deple- 
tion-accumulation region where u, < u H is given in Fig. 31(a). The 
equivalent circuit in the mid-gap region where u, = u B ± o few kT/q 
is shown in Fig. 31(b). Finally, the equivalent circuit in the weak 
inversion region where u B < u. < 2u n is shown in Fig. 31(c). 

(x) It is impossible to make an MIS structure in which surface po- 
tential is absolutely uniform because even if built-in charge and inter- 
face state density is made small, random distribution of ionized bulk 
impurities or random fluctuations of oxide thickness will still be present. 

(xi) The accuracy required for extracting this detailed a picture of 
the Si-Si0 2 interface from MIS capacitor measurements could only be 
provided by the conductance technique (see Fig. 1). 



■ 




(c) 



Fig. 31 — (a) Schematic of interface state branch of equivalent circuit in deple- 
tion-accumulation region where u, < ub. (b) Schematic of interface state branch 
of equivalent circuit in mid-gap region where u. = u B ± a few kT/q. (c) Schematic 
of interface state branch of equivalent circuit in weak inversion region where 
ub < u, < 2ub . 
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APPENDIX A 

In this Appendix, interface state admittance is derived from (3) : 

i.(t) = qN.c n [l - Kt)]n.(t) - qN,e n j(t)- (65) 

Expressing f(t) as the sum of a dc and an ac part: 

m = /. + */> ( 66 > 

where /„ is the Fermi function established by the bias and 8f the change 
caused by the ac signal. Similarly, 

n.{t) = n. + 8n, , (67) 

where n,„ is the electron density at the silicon surface when the Fermi 
level is at the trap level and 8n, the change caused by the ac signal. 
Substituting (66) and (67) into (65) and making the small signal ap- 
proximation by neglecting second order terms, we get 

i t (t) = qN,c n [(l - />.„ + (1 - jo) 5n. ~ n t0 8f] - qe n N,(1„ + 5/). 

(68) 

From detailed balance 

qN.c n tt ~ />.. = Q^N.1„ . (69) 

Solving (69) for e„ and substituting into (68) we get 



(1 - /.) 8n. - n lo | j 



i,(t) = qN,c n 
The net current density can also be expressed as 



(70) 



m = qN.§. (71) 
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Equating (71) to (70) 

% = c(l -/„)5n. -c„n..f- (72) 

at U 

The small signal variation of the Fermi function is 5/ = / m exp (jut). 
Combining this with (66) results in 

I = * * ™ 

Substituting (73) into (72) and solving for 5/ 

1 " n..(l + jaf /c n n..) K J 

Substituting (74) into (70) results in 

l,(() -(iT^M.)n.; (75) 

The ratio Sn,/n l0 is equal to the ac surface potential. This can be 
seen as follows: 

8n, w, exp (m. — M fl ) 8u, _ ,-gx 

n IO " w, exp (w. — ii B ) 

where n, is the intrinsic carrier density and u, is the silicon band bending 
or surface potential in units of kT/q and u B the potential difference 
between mid-gap and Fermi level in the quasi neutral region of the 
silicon in units of kT/q. Expressing 8u, in volts is 

«u. - ^ «*. . (77) 



Equation (75) becomes 



t.(0 = J" fcr (1 + Ja)/o/Cnr0 • 



(78) 



This can be written as 

t,(0 = 7. 5-A. , (79) 

where 

^• = p> kf{l+Mc n n. ) (80) 
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APPENDIX B 

B.l Random Fluctuations oj Impurity Charges in the Silicon 

We shall use the model proposed by Shockley 23 for this calculation. 
Otherwise the procedure for finding the standard deviation of the 
surface potential is similar to that for the case of the surface charges. 
Let M be the number of ionized acceptor* in a characteristic volume 
in the silicon space-charge region. When M is large, the probability 
that there are M ionized acceptors in a characteristic volume is given 
by the Gaussian approximation of a Poisson distribution 

P(M) = (2ttM)- 5 exp [-(M - Mf/2M]. (81) 

The characteristic volume from Ref. 23 is a cube in which the length 
of a side is equal to the space-charge width W. The first transformation 
of the probability density function is from number of charges to number 
density 

P(N A ) = p(M) dM/dN A , (82) 

where N A is the ionized acceptor density. The relation between M 
and N A is 

M = N A W a . (83) 

Combining (81), (82), and (83), the result of the first transformation is 

P(N A ) = (2wN A WT h W' A exp [-W\N A - N A )*/2N A l (84) 

where N A is the mean acceptor density. 

The second transformation from number density to surface potential is 

P(u.) = P(N A ) dNJdu, . (85) 

Rewriting (31) and (32) results in 

Q = C ox (v + UM + (2qe Bi N A /^[exp (-u.) + U. - 1]K (86) 

Solving (86) for N A is 

Na ~ \2qeJ [exp (-w.) + W. - 1] K ° } 

Equation (87) is not a single-valued relation between N A and u t . We 
therefore restrict the problemas with the surface charges to the case 
where the fluctuations N A — N A are very small. For very small fluctua- 
tions, (87) can be differentiated assuming Q and oxide thickness to be 
uniform to get 
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dN = J(Ma)* c„ 

Wqt.iP/ [exp (-m,) + u. - 1] J 

iV,[l-exp(-u,)1 | 
|exp (— u,) + w. — 1]J v ' 

Equation (88) is now expressed as a function of u, using (32) 

c - » - ' Sr - [**i5SfenJi> - - <-•»• « 

Solving (89) for exp (— w.) + w. — 1, we have 

exp (-u,) + it. - 1 ■= ^^ ]^ 2 [1 - exp (-w.)] 2 . (90) 

Substituting (90) into (88) we get 

<//31r [1 — exp (—«,)] v ' 

Because A' x is given by (87) when u, = u, , dN A can be evaluated 
about A r A at each bias from 

,, r 2[W(u,)C ox + e.,1 

dN * = — Wifk^n i - ^du, . (92) 

qpw (w,)[l — exp (— w,)] v ' 

Equation (92) is the transformation equation sought. Replacing dN A 
and du, in (92) by the small fluctuations N A — K' A and u, — u, re- 
spectively, we get 

*r tf 2[W(tZ I )C„ + e./l 

<//3lr (m,)[1 — exp (— -u,)] v ' K ' 

Combining (84), (85), (92), and (93), we get 

P(u.) = (2**1)-'' exp [-(«. - u.Y/2<r%l (94) 

where the standard deviation a„ is 

q(3[ N A W(u.)ftl ~ exp (-U.)] 

2[W(u,)C oz + «„•] " ( 95 > 

B.2 Random Fluctuations of Oxide Thickness 
Solving (80) for x, the oxide thickness, 

is. « « . Up. + «./0) , Qr . 

C„ Q - ((2g«.,JV y( //3)[exp (-«.) + u . - 1]}* ' W 
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where e„ is the dielectric permittivity of the oxide layer. If Q is not 
zero, it will influence the fluctuations of surface potential caused by 
fluctuations in oxide thickness even if it is uniform. To get the influence 
of oxide thickness fluctuations independent of surface charge, we let 
Q = 0. Thus, (96) becomes 

_ e > o + U./0) (q7) 

X ~ K2ff€.^//5)[exp(-«.)+u. -1])* V ' 

Differentiating (97), assuming N A to be uniform, results in 



dx (2qe. i N A /^ Wexp (-u.) + u. - l] 1 

_ (v + u./m- e*P (-u.)] \ d m 

2[exp (-w.) + m. - 1]' J 

Eliminating exp (— «.) + u. - 1 from (98) by using (90) we get 

dx = — a , T w2ri -, ^ du, . (99) 

q(3N A W [1 — exp (-w.)J 

Equation (97) gives x when u, = u. so that dx evaluated about x is 

W(u.)e ox + fe,, ,_ /mm 

dx is the standard deviation of x and dw, the standard deviation of 
surface potential thus, 

_ = ^^.)[l - exp (-«.)] ^ 
Llr (M,)e 0I + xe.,J 

APPENDIX C 

There is a simple way to relate C,(w, m.) to measured (GJu>) mn . 
The distributed interface state network in Fig. 8(a) can be represented 
at one frequency and bias by an equivalent series R.C, network in 
which R, and C. are frequency and bias dependent. The equivalent 
circuit will then look like Fig. 6. If we let R, and C, be independent 
of frequency, G p /oj vs frequency from (17) can be made to coincide 
with the measured GJu vs frequency curve at each bias at (G P /o)) maXi 
as shown in Fig. 32. Therefore, from (17), C.(oj, u.) = 2(G p /u)) mux at 
each bias and we obtain (54) by substituting this into (53) and solving 
for C D . 



ELECTRICAL PROPERTIES OF SILICON-OXIDE INTERFACE 



1127 




LOG FREQUENCY 



Fig. 32 — Schematic of G,/io vs log frequency showing how C, is obtained 
from measured G P /w peak. 



APPENDIX D 



We shall calculate G p /co vs oj which would result from transitions 
between the conduction band near the surface and states distributed 
with uniform density into the oxide. The time constant distribution 
for tunnelling from (61) is 



r, = t 0< exp {2K&, 



(102) 



where t„, = (n, va ) \ 

The contribution to GJu from states located between £ and £ + tf£ 
in the oxide is 

d(GM = (q/2u>)(N../x L )T. l exp (-2K.Q 

• In [1 + u 2 tI, exp (4K£)] d$, (103) 

where (102) has been substituted for t in (20) for the continuum of 
states, N„/x L is the uniform density of states in the oxide, cm~ 3 -eV~ l , 
and x L is the oxide thickness, cm. 

Total GJu from all the states in the oxide therefore, is 

G> = (q/2»XN../x L )T., r exp (-2K£ 

Jo 

•In [1 + w 2 r 2 0< exp (4K£\ d£. (104) 
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To get GJw vs w, (104) is numerically integrated on an IBM 7094 
computer using the trapezoidal rule. The result is shown in Fig. 33. 
Because we are interested only in the shape of the curve, the frequency 
scale is arbitrary. 

GJo) is normalized to eliminate N„/x L . A plot of normalized G p /u 
obtained from experiment is also shown in Fig. 33 for comparison. 

The shape of the curve calculated from the tunnelling model is 
determined essentially by (102). The shortest time constant in this 
distribution is when £ = which is right at the interface. GJw drops 
off rapidly as the period of the applied signal becomes comparable 
to and shorter than this time constant. When the period of the applied 
signal is long compared to the shortest time constant in (102), G,/u 
becomes a constant as we have assumed a constant density of states 
into the oxide. The main purpose of assuming a uniform density of 
states into the oxide is to illustrate the asymmetry of the GJu vs log a> 
curve introduced by the tunnelling time constant distribution (102). 
The GJw vs log a> curve will go through a maximum if the density of 
states decreases with distance into the oxide. As long as the density 
of states decreases more slowly with distance into the oxide than the 
time constants given by (102) increase, GJo will decrease faster on 
the high frequency side of the peak than on the low-frequency side. 




FREQUENCY (ARBITRARY SCALE) 



Fig. 33 — Normalized G P /u vs log frequency for tunnelling calculated from 
(104) and normalized G P /w vs log frequency obtained for measurement. Frequency 
scale is arbitrary to compare the shapes of two curves. 
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The G„/u vs log w curve obtained from measurement is symmetric 
about its peak. A distribution of states into the oxide which would 
would give a G p /u vs log u curve symmetric about its peak using the 
tunnelling model would be a peaked density distribution having its 
maximum some distance into the oxide. In addition, the width ofthis 
distribution would have to be considerably smaller than 1/2K to get 
G p /u) vs log a) curves in the frequency range between 50 Hz and 500 kHz. 
Such a distribution would be rather arbitrary. Thus, it can be seen 
from the symmetry of the observed G p /w vs log o curves that tunnelling 
to states distributed into the oxide is not likely to be the dominant 
cause of time constant dispersion in the samples measured. 

LIST OF SYMBOLS 

a Maximum amplitude of the ac surface potential in volts. 

b Maximum amplitude of applied ac signal in volts. 

C Equivalent parallel capacitance of interface state branch 

of equivalent circuit having a single time constant in 

farads/cm 2 . 
CB Representation of conduction band. 

C D Depletion layer capacitance in farads/cm 2 . 

C t Inversion layer capacitance in farads/cm 2 . 

C m Measured equivalent parallel capacitance of the MIS 

capacitor in farads/ cm 2 . 
C Electron capture probability in cm 3 /sec. 

C om Low-frequency capacitance of MIS capacitor in farads/cm 2 . 

C ox Oxide layer capacitanec in farads/cm 2 . 

c p Hole capture probability in cm 3 /sec. 

C„ Equivalent parallel capacitance of depletion layer and 

interface state branch of equivalent circuit in farads/cm 2 . 
C. Interface state capacitance associated with a single time 

constant in farad/cm". 
dv Number of characteristic areas in which the number of 

surface charges is between N and N -f- dN or the surface 

potential is between u, and u, + da, . 
E F Fermi level in units of IcT/q. 

E { Mid-gap energy in units of kT/q. 

e n , e p Emission constants for electrons and holes in sec" . 

/„, Maximum value of the perturbation of the Fermi function 

by the ac signal. 
/„ Fermi function at the dc bias. 
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f(t) Fermi function as a function of time. 

G m Measured equivalent parallel conductance of MIS capac- 

itor in mhos/ cm 2 . 

G n {t) Emission rate of electrons from a single level interface 

state in cm~ z -sec~ l . 

G v Equivalent parallel conductance of the interface state 

branch of the equivalent circuit in mhos/cm 2 . 

G' v Equivalent parallel conductance corrected for C ox in mhos. 

O p (t) Emission rate of holes from a single level interface state 

in cm" 2 -sec -1 . 

G. Majority carrier capture conductance of a single level 

interface state in mhos/cm 2 . 

i, Current density charging the interface states in amp/cm 2 . 

i, c Current density charging the space-charge layer in amp/cm . 

i T Current density through MIS capacitor in amp/cm 2 . 

h Boltzman's constant in eV-coul/°K. 

M Number of ionized acceptors in a characteristic volume. 

M The mean of M . 

N Number of built-in charges and charged interface states 

in a characteristic area. 

N The mean of N. 

N A Acceptor impurity density in cm 

N A Mean ionized acceptor density in cm 

N D Donor impurity density in cm 

rii Intrinsic carrier concentration in cm 

N. Density of a single level interface state in cm 

n,„ Electron density established by the dc bias in cm -3 . 

N„ Density of interface states, dN,/d\p in cm" 2 -eV _1 . 

n,{t) Time dependent electron density at the silicon surface 



-3 

i en 

-■'■ 

-3 

-2 



in cm 3 . 



P(M) Probability that there are M ionized acceptors in a char- 

acteristic volume. 

P(N) Probability of finding N charges in a characteristic volume. 

P(Q) Probability that the surface-charge density is Q in cm 2 /coul. 

p.„ Hole density established by the dc bias in cm -3 . 

p,(t) Time dependent hole density at the silicon surface in cm 3 . 

P(u,) Probability that surface potential in a characteristic area 

is u, . 

q Electronic charge in coulomb. 

Q Sum of built-in and interface state charge densities in 

coul/cm 2 . 
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Q Mean surface-charge density in coul/cm 2 . 

Q, Fixed charge density in the oxide in coul/cm . 

Q, Total interface state charge density in coul/cm 2 . 

Q, e Silicon space-charge density in coul/cm 2 . 

Q T Total charge density in coul/cm 2 . 

R„_, ; R Pi , Electron and hole capture resistance of a single level in 
ohms-cm 2 . 

R„(t), R p (i) Capture rate of electrons and holes by a single level interface 
state in cm~ 2 -sec _1 . 

R, Majority carrier capture resistance of a single level inter- 

face state in ohm-cm 2 . 

t time in seconds. 

T Absolute temperature in °K. 

u Potential difference between mid-gap and an interface 

state level in units of kT/q. 
= In (rii/Nn), for n-type silicon 



{= In (N A /ni), for p-type silicon. 

v P Potential difference between mid-gap and Fermi level at 

silicon surface in units of kT/q 

u. In w/c„n,- . 

u, Surface potential in units of kT/q. 

u. Mean surface potential in units of kT/q. 

VB Representation of valence band. 

v Average thermal velocity of holes and electrons in cm/ sec. 

v a (t) Time dependent voltage applied to MIS capacitor in volts. 

v. dc bias across MIS capacitor in volts. 

v. x Arbitrary dc bias in volts. 

W Depletion layer thickness in cm. 

x Mean oxide thickness measured from the interface in cm. 

?/ In wr m = u B + u — u. . 

Vm u, — u. + In 2.5. 

Y, Admittance of a single level interface state in mhos/cm . 

Y„ Admittance of an interface state continuum in mhos/cm 2 . 

Z, MIS capacitor impedance less reactance of the oxide layer 

in ohms-cm 2 . 

a Characteristic area, cm 2 

|8 = 38.5 q/kT in volt -1 . 

A Additive constant in volts. 

5/ Time dependent part of Fermi function. 

bn. Incremental electron surface-charge density induced by 
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5v a ac signal applied to MIS capacitor in volts. 

8\f/ M Silicon band bending by ac perturbation in volts. 

e fl , e ux Dielectric permittivity of silicon and Si0 2 layer respectively 
in farad/cm. 

t Distance into oxide measured from interface in cm. 

a B Standard deviation of surface potential caused by fluctua- 

tions of ionized acceptors. 

<r n , a p Electron and hole capture cross sections in cm 2 . 

a. Standard deviation of surface potential caused by fluctua- 

tions of surface charges. 

a z Standard deviation of surface potential caused by oxide 

thickness fluctuations. 

r Majority carrier time constant of a single level interface 

state, R,C, in seconds. 

Tm Measured interface state time constant in seconds. 

T mo Time constant at u. = u B . 

t„,,;t p ,. Electron and hole time constant for a single level in seconds. 

\f/ Energy in electron volts. 

\J/ g (kT/q) u B in volts. 

y]/. Band bending or surface potential in volts. 

$, Mean surface potential in volts. 

\f/, t Surface potential plus an additive constant in volts. 

}J/ to Surface potential at the dc bias in volts. 

a; Angular frequency = 2irj in sec . 
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